﻿ Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Conducător ştiinţific: Prof dr ing Miranda Naforniţă Doctorand: as ing Horia Balta Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta Cuprins 1 Sisteme de transmisiuni cu spectru împrăştiat 1 1 1 Principiul împrăştierii spectrului Modelul unui sistem de comunicaţie cu spectru împrăştiat 1 1 2 Avantajele comunicaţiilor cu spectru împrăştiat Aplicaţii 3 1 3 Tehnici de împrăştiere a spectrului Comparaţie 12 1 3 1 Împrăştiere cu secvenţa directă 12 1 3 2 Împrăştiere cu salt de frecvenţă 14 1 3 3 Împrăştiere cu salt de timp 14 1 3 4 Sisteme hibride 15 2 Secvenţe pseudo-aleatoare (CN—coduri zgomot) 17 2 1 Cerinţele codurilor zgomot 17 2 2 Secvenţe pseudo-aleatoare Proprietăţi Generare 18 2 2 1 Secvenţe Registru de Deplasare 18 2 2 2 Proprietăţile secvenţelor M 21 2 3 Coduri utilizate pentru împrăştierea spectrului 23 2 3 1 Coduri Walsh Hadamard 23 2 3 2 Secvenţe M 26 2 3 3 Secvenţe Gold şi Kasami 26 2 3 4 Soluţii alternative Concluzii 28 3 Tehnici adiţionale de îmbunătăţire a performanţelor sistemelor cu spectru împrăştiat 30 3 1 Întreţeserea 30 3 2 Codarea corectoare de erori 32 3 2 1 Decodare soft / hard 32 3 2 2 Câştigul de codare 32 3 2 3 Criteriul MAP 33 3 3 Concatenarea 34 3 3 1 Concatenarea codurilor 34 3 3 3 Coduri bloc produs 34 3 3 4 Coduri convoluţionnale concatenate 35 4 Coduri grup 38 4 1 Coduri grup 38 4 2 Coduri Reed–Muller (R-M) 39 5 Coduri ciclice 43 5 1 Coduri ciclice –descriere generală 43 5 2 Codarea codurilor ciclice 44 5 2 1 Registru de Deplasare cu Reacţie 44 5 2 2 Codor ciclic cu RDR şi sumatoare exterioare 46 5 2 3 Codor ciclic cu sumatoare multiple 47 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta 5 3 Decodarea codurilor ciclice 49 5 3 1 Decodor cu RDR pentru un cod ciclic detector de erori 50 5 3 2 Decodor cu RDR pentru un cod ciclic corector de o eroare 54 5 4 Coduri BCH 59 5 4 1 Construcţia cuvintelor de cod BCH 59 5 4 2 Decodarea codurilor BCH 60 5 4 3 Codul Golay 61 5 5 Coduri Reed–Solomon (R-S) 62 q 5 5 1 Câmpul Galois GF(2) 62 5 5 2 Polinomul generator al codului 64 5 5 3 Decodarea codurilor R-S 65 6 Coduri convoluţionale 69 6 1 Coduri convoluţionale –descriere generală 69 6 1 1 Codor convoluţional 69 6 1 2 Diagrama trellis 70 6 1 3 Distanţa de cod 71 6 2 Decodor maximum plauzibil –algoritmul Viterbi 73 6 2 1 Decodarea MAP 73 6 2 2 Metrica ramurii 73 6 2 3 Algoritmul Viterbi 75 6 3 Generalizarea exemplului precedent 77 Bibliografie 79 Anexa A Lista polinoamelor primitive Anexa B Perechi de polinoame ce generează secvenţe preferate Anexa C Corespondenţa resturi –cuvinte eroare corectabile Anexa D Lista polinoamelor ireductibile Anexa E Program de simulare a unei scheme de HCCC Anexa F Decodarea în frecvenţă a codurilor R-S Algoritmul Berlekamp–Massey Anexa G Exemplificarea algoritmului Viterbi Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta 1 Sisteme de transmisiuni cu spectru împrăştiat 1 1 Principiul împrăştierii spectrului Modelul unui sistem de comunicaţie cu spectru împrăştiat Împrăştierea spectrului este o tehnică de modulaţie a semnalului purtător de informaţie Această tehnică reprezintă o soluţie calitativ superioară în multe aplicaţii ce implică comunicaţii prin medii cu acces multiplu, în special canalele radio Spre deosebire de tehnicile clasice de modulaţie, împrăştierea spectrului aduce câteva elemente noi Primul este că semnalele cu spectru împrăştiat utilizate pentru transmiterea informaţiei digitale se disting prin aceea că lăţimea lor de bandă W [Hz] este mult mai mare decât rata informaţiei R [biţi per secundă] Astfel factorul expansiunii de bandă β e = W/R (1 1) pentru un semnal cu spectru împrăştiat este mult mai mare decât unitatea Redundanţa mare, inerentă în semnale cu spectru împrăştiat, este necesară pentru a putea realiza transmisii la nivelele mari de interferenţă ce sunt întâlnite în unele canale radio sau pe satelit Codarea este un element important în proiectarea semnalelor cu spectru împrăştiat deoarece formele de undă codate sunt de asemenea caracterizate printr-un factor de expansiune a benzii ce este mult mai mare decât unitatea şi deoarece codarea este o metodă eficientă de introducere a redundanţei Al doilea element important implicat în proiectarea semnalelor cu spectru împrăştiat este pseudoaleatorul, care face ca semnalele să apară similar cu zgomotul aleator şi dificil de demodulat de receptoare, altele decât cel dedicat Acest element este descris împreună cu aplicaţia acestor semnale Semnalele cu spectru împrăştiat sunt utilizate cu scopul: (1) de a combate sau suprima efectele dăunătoare ale interferenţei datorate bruiajului, a interferenţei cu alţi utilizatori ai canalului şi a interferenţei proprii datorate propagării pe căi multiple; (2) de a ascunde semnalul prin transmiterea sa la o putere joasă şi, astfel, de a face dificilă detectarea sa de către un ascultător neavizat în prezenţa zgomotului de fond; (3) crearea de mesaje private în prezenţa altor ascultători (4) În alte aplicaţii decât comunicaţiile, semnalele cu spectru împrăştiat sunt folosite pentru a obţine un nivel de acurateţe (timp de întârziere) cum şi un nivel al ratei (vitezei) cum măsurătorilor în radar sau navigaţie În combaterea bruiajului, este important pentru utilizatori că bruiorul, care încearcă să perturbe comunicaţiile, nu cunoaşte dinainte caracteristicile semnalului, cu excepţia întregii benzi a canalului şi a tipului modulaţiei (PSK, FSK, etc) care sunt folosite Dacă informaţia digitală este codată „clasic”, un bruior sofisticat poate imita uşor semnalul emis de transmiţător şi, aşadar, să deruteze receptorul Pentru a preveni această posibilitate, 1 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta transmiţătorul introduce un element de impredictibilitate sau aleator (pseudoaleator) în fiecare formă de undă a semnalului transmis codat care este cunoscut de receptorul avizat dar nu şi de bruior Ca o consecinţă, bruiorul trebuie să sintetizeze şi să transmită un semnal de interferenţă fără a cunoaşte forma pseudoaleatoare Interferenţa cu alţi utilizatori apare în sistemele de comunicaţie cu acces multiplu în care un număr de utilizatori împart o bandă comună a canalului La orice moment dat, un subset al acestor utilizatori transmit simultan informaţie prin canalul comun la receptorii corespunzători Presupunând că toţi utilizatorii folosesc acelaşi cod pentru codare şi decodare pentru secvenţele lor de informaţie, semnalele transmise în acest spectru comun pot fi distinse 1 unul de altul prin impunerea unei forme pseudo-aleatoare, de asemenea numită cod, în fiecare semnal transmis Astfel un receptor particular poate descoperi informaţia transmisă şi adresată lui prin cunoaşterea formei pseudoaleatoare, altfel spus a cheii, utilizate de transmiţătorul corespondent Acest tip de tehnică de comunicaţie, care permite utilizatorilor multipli să folosească simultan un canal comun pentru transmiterea informaţiei, se numeşte acces multiplu cu diviziune în cod (CDMA—Code Division Multiple Access) Componentele reductibile multicale, rezultate din propagarea dispersivă în timp prin canal, pot fi văzute ca o formă de interferenţă proprie Acest tip de interferenţă poate fi de asemenea suprimată prin introducerea formei pseudoaleatoare în semnalul transmis, aşa cum va fi descrisă ulterior Un mesaj poate fi ascuns în zgomotul de fond prin împrăştierea benzii sale prin codarea şi transmiterea semnalului rezultat la o putere mică Datorită nivelului scăzut al puterii sale, semnalul transmis se spune că este ”ascuns” (“covert”) Are o probabilitate scăzută de a fi interceptat (detectat) de un ascultător ne-avizat, şi, de aceea, este numit semnal cu probabilitate mică de intercepţie (LPI) În final, mesajul privat poate fi obţinut prin impunerea unei forme pseudo-aleatoare în mesajul transmis Mesajul poate fi demodulat de către receptorii avizaţi, ce cunosc forma pseudo-aleatoare sau cheia folosită la transmiţător, dar nu şi de alţi receptori ce nu au cunoştinţă despre cheie Schema bloc prezentată în Figura 1 1 ilustrează elementele de bază ale unui sistem de comunicaţie digital cu spectru împrăştiat, cu o secvenţă de informaţie binară la intrarea transmiţătorului şi la ieşirea receptorului Codorul şi decodorul de canal, cât şi modulatorul şi demodulatorul sunt elementele de bază ale sistemului În completarea acestor elemente avem două generatoare identice de secvenţe (tipar) pseudoaleatoare, una ce conlucrează cu modulatorul la finele transmisiei iar cea de-a doua ce conlucrează cu demodulatorul la finele recepţiei Generatoarele generează o secvenţă binară pseudoaleatoare sau pseudonoise (PN) ce este introdusă în semnalul transmis la modulator şi extrasă din semnalul recepţionat la demodulator Sincronizarea secvenţei PN generate la receptor cu secvenţa PN conţinută în semnalul recepţionat este necesară pentru demodularea semnalul recepţionat Iniţial, înainte de transmiterea informaţiei, sincronizarea poate fi obţinută prin transmiterea unui tipar de bit pseudoaleator pe care receptorul îl va recunoaşte, şi în prezenţa interferenţei, cu mare probabilitate După realizarea sincronizării generatoarelor, poate să înceapă transmisia informaţiei Interferenţa este introdusă în transmisie la trecerea semnalului de informaţie prin canal Caracteristicile semnalului ce interferă depind în mare măsură de originea sa Interferenţa poate fi clasificată ca fiind fie de bandă largă fie de bandă îngustă relativ la lăţimea de bandă a semnalului purtător de informaţie; de asemenea, semnalul de interferenţă poate fi continuu în timp sau discontinuu (pulsatoriu) în timp De exemplu, un semnal de 1 în literatură sunt întâlnite denumirile: cod zgomot sau pseudo-zgomot (PN—pseudo noise), semnătură sau cheie În lucrarea de faţă vor fi utilizate, după caz, toate aceste denumiri 2 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta bruiaj poate consta dintr- una sau mai multe sinusoide în banda de transmisie utilizată pentru transmiterea informaţiei Frecvenţele sinusoidelor pot rămâne fixe sau se pot schimba cu timpul Ca un al doilea exemplu, interferenţa generată în CDMA de alţi utilizatori ai canalului poate fi de asemenea de bandă largă sau îngustă, depinzând de tipul semnalului cu spectru împrăştiat utilizat la realizarea accesului multiplu Dacă este de bandă largă, poate fi caracterizată ca un zgomot echivalent, gaussian, alb, aditiv Secvenţa de Datele de informaţieş ie ire Codor de ModulatorDemodulatorDecodor Canal de canal canal Generator de Generator de tipar tipar pseudo-aleatorpseudo-aleator Figura 1 1 Model de sistem digital de comunicaţie cu spectru împrăştiat Modulaţiile utilizate în sistemele de transmisiuni cu spectru împrăştiat pot fi de tip PSK sau FSK Modulaţia PSK este adecvată aplicaţiilor unde coerenţa de fază dintre semnalul transmis şi semnalul recepţionat poate fi menţinută într-un interval de timp ce este relativ lung comparativ cu 1/W, unde W reprezintă banda semnalului transmis Pe de altă parte, modulaţia FSK este potrivită aplicaţiilor în care o astfel de coerenţă de fază nu poate fi menţinută, datorită variaţiei în timp a parametrilor semnalelor din canalele de comunicaţii Acesta poate fi cazul unei legături de comunicaţie dintre două avioane de mare viteză sau dintre un avion rapid şi un terminal de pe sol Secvenţa PN generată la modulator este utilizată împreună cu modulaţia PSK pentru a deplasa faza semnalului PSK pseudo-aleator Semnalul modulat rezultat este numit semnal cu spectru împrăştiat cu secvenţă directă (DS) sau pseudonoise (PN) Utilizată împreună cu FSK binar sau M-ar, secvenţa pseudoaleatoare selectează frecvenţa semnalului transmis în mod pseudoaleator Semnalul rezultat este denumit semnal cu spectru împrăştiat cu salt de frecvenţă (FH) 1 2 Avantajele comunicaţiilor cu spectru împrăştiat Aplicaţii În acest paragraf vor fi prezentate câteva dintre aplicaţiile ce utilizează semnale cu spectru împrăştiat a)-Estimarea întârzierii semnalului reflectat în sistemele radar Dispersia timpului de întârziere are expresia aproximativă: 2 21τ c σ =≈, (1 2) 2222τ ()() 2W2πξπξ 3 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta 2 unde: ξ = raportul semnal pe zgomot; W = banda semnalului; τ = 1/W = dublul timpul de c creştere Relaţia (1 2) arată că, pentru a avea o bună precizie în estimare, este necesar să se utilizeze un semnal de bandă largă Dacă durata semnalului, T, este aleasă mică, atunci este 2 necesar (în vederea unui raport semnal pe zgomot ξ mare) un maxim de putere O alternativă este alegerea lui T>>τ c, sau echivalent: W T >> 1, (1 3) relaţie satisfăcută prin utilizarea semnalului cu spectru împrăştiat Raţionamentul şi concluzia sunt identice cu cele prezentate în cazul în care se doreşte o estimare a frecvenţei semnalului b)-Obţinerea de răspunsuri scurte şi cu RSZ mare la ieşirea filtrelor adaptate (FA) sau corelatoarelor Filtrul adaptat poate separa în timp semnalele de intrare, ce sunt întârziate şi suprapuse unele peste altele şi cu zgomotul de fond, dacă timpul de corelare, τ c, este mai mic decât distanţa în timp dintre ele, Figura 1 2 Condiţia ca τ Ti este realizată prin c > 1 g)-Combaterea fading-ului Fading-ul (fluctuaţia semnalului recepţionat) este un efect al propagării semnalului pe mai multe căi, de la emiţător spre receptor, în cazul transmisiilor radio terestre Propagarea diferitelor părţi din semnalul emis de antena emiţătoare spre antena receptoare se face, pe lângă „unda directă”, prin reflexii multiple Astfel, rezultatul este că la recepţie sosesc diferite semnale, copii întârziate şi/sau atenuate ale undei directe, suprapuse peste aceasta În plus aceste întârzieri şi atenuări variază în timp Urmarea este că, la aceeaşi putere a semnalului emis, probabilitatea de eroare la reconstrucţia semnalului în receptor creşte dramatic dacă componentele multicăi se suprapun în antifază, Combaterea fading- ului se realizează prin „diversitate” Principial, diversitatea presupune existenţa unui număr de L canale paralele, purtând aceleaşi date, pentru care fading-ul este statistic independent rii complete a semnalelor din toate Datorită independenţei probabilităţii, PL, a suprimă L canalele faţă de a unuia, aceasta va descreşte cu L: PL = P, unde P0 este probabilitatea 0 suprimării semnalului într-un canal Modurile „clasice” de obţinere a diversităţii sunt: –diversitate spaţială, prin dispunerea mai multor antene de recepţie, separate între ele prin câteva lungimi de undă Fading-urile în ele sunt independente –diversitate în frecvenţă, prin transmisia paralelă a mai multor frecvenţe purtătoare cu diferite lungimi de undă Ambele metode fac uz de diferenţa de fază dintre diferitele căi: ∆ϕ = 2⋅π⋅ ∆D / λ, (1 9) 6 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta D este diferenţa de lungime dintre căi iar λ este lungimea de undă Astfel, considerând unde ∆ semnalele reprezentate prin vectori, este posibil ca o antenă de recepţie (vorbind despre cazul diversităţii spaţiale) să compună vectorii—semnal—recepţionat ca şi în Figura 1 5 a, (cazul unei suprimări a semnalului), dar este mult mai probabil ca să existe măcar o antenă receptoare (atunci când sunt mai multe) care să compună vectorii—semnal—recepţionat ca şi în Figura 1 5 b, caz în care avem o amplificare a semnalului şi deci o recepţie posibilă Asemenea, şi în cazul diversităţii în frecvenţă, pentru unele lungimi de undă este posibil să existe suprimare (Figura 1 5 a), dar este mult probabil ca măcar pentru o lungime de undă să existe amplificare Si S S = Σ 1 Sn S1S2S3 Si S2 S = Σ Sn S3 a) suprimare; b) amplificare; Figura 1 5 Diagrama fazorială în cazul recepţiei prin „diversitate” O alternativă la metodele prezentate anterior o reprezintă diversitatea temporală Aceasta se bazează pe separarea temporală a semnalelor diferitelor căi, la recepţie, prin filtrare adaptată Diferenţa de timp dintre semnalele aferente a două căi este ∆τ = ∆D / c, (1 10) a lui ∆τ peste toate unde c este viteza luminii Astfel, dacă notăm cu ∆τ min valoarea minimă spunsurile filtrului perechile de căi, şi dacă ∆τ min > 2/W = 2⋅τc, atunci conform Figurii 1 2, ră adaptat la semnalele aferente diferitelor căi pot fi separate în timp După aceasta, ele pot fi procesate într-o manieră coerentă sau necoerentă Dacă fiecare amplitudine complexă, Ai, poate fi măsurată suficient de precis, atunci cea mai bună soluţie este să sumăm semnalele, * ponderate cu Ai (procesarea coerentă ): L ∗ () () ts, ⋅=tsA (1 11) ∑ iii 1i = ii i În acest fel, în schimbul pierderilor, vom avea o creştere unde si(t) este semnalul aferent că 2 în RSZ de L ori, iar utilizarea spectrului împrăş tiat (cea care permite separarea căilor) face ca efectul propagării multicale să fie unul benefic Combaterea fading-ului se poate face combinând diversitatea cu alte metode, cum ar fi codarea, 7 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta h)-Accesul multiplu Problema accesului multiplu apare în situaţiile în care un număr relativ mare de utilizatori, poziţionaţi într-o zonă geografică comună, utilizează, (posibil) 1 simultan, aceeaşi alocaţie spectrală în vederea comunicării separate sau private În abordarea clasică, fiecare utilizator al sistemului cu acces multiplu este înzestrat cu resurse certe, ca frecvenţa (în tehnica numită FDMA –Frequency Division Multiple Access) sau sloturile de timp (în tehnica numită TDMA –Time Division Multiple Access), sau ambele, care sunt disjuncte cu cele ale oricărui alt utilizator În acest fel canalul cu acces multiplu se reduce la o multititudine de canale punct la punct, asigurând o izolare perfectă a resurselor fiecărui utilizator de a celorlalţi S1(ƒ)s1(t) WuTu ƒ t S2(ƒ)s2(t) t ƒ W=M⋅WuT=M⋅ Tu SM(ƒ)sM(t) t ƒ TDMA–fiecare utilizator ocupă FDMA –fiecare utilizator ocupă tot toată banda disponibilă şi doar o timpul disponibil şi doar o parte 1/M parte (slot) 1/M din timp; din banda disponibilă; Figura 1 6 Sisteme de comunicaţie cu acces multiplu „tradiţionale” Tu TDMA FDMA WW M M 1 2 Wu 2 Resursa de bandă Resursa de bandă 1 de frecvenţă de frecvenţă Resursa de timp TResursa de timp T Tu≥1/W⇒M≤WTWu≥1/T⇒M≤WT Figura 1 7 Alocarea resurselor şi numărul de utilizatori 1 Apelativul „separate” subliniază faptul că nu este vorba de o comunicaţie prin „difuzare”, caz în care un emiţător transmite către mai mulţi utilizatori; comunicaţiile invocate în text reprezintă legături multiple între diverşi emiţători şi receptori 8 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta Capacitatea fiecărui canal este limitată numai de lăţimea de bandă şi timpul alocat, de degradarea cauzată de zgomotul de fond, în principal de origine termică, şi de anomaliile de propagare, care produc efectele de fading multicăi (multipath fading) şi de umbrire (shadowing effect) Astfel, în aparenţă, capacitatea, compusă de componentele individuale ale canalelor cu acces multiplu, va fi egală cu capacitatea unui singur utilizator ocupând întregile resurse compuse ale canalelor individuale O alternativă a metodelor clasice, prezentate anterior, este CDMA—Code Division Multiple Access Această metodă alocă (vezi Figura 1 8) toate resursele simultan tuturor utilizatorilor, iar puterea transmisă de fiecare este cea minim necesară pentru a menţine un raport semnal pe zgomot dat de nivelul de performanţă necesar Fiecare utilizator întrebuinţează un semnal cu lăţime de bandă ca cea a zgomotului ocupând întreaga bandă alocată, pentru cât timp e necesar În acest fel fiecare utilizator contribuie la zgomotul de fond ce afectează toţi utilizatorii, dar în modul cel mai puţin distructiv Această interferenţă adiţională limitează capacitatea, dar pentru că resursele alocate de timp şi lărgime de bandă sunt nerestricţionate, capacitatea rezultantă este semnificativ mai mare decât la sistemele convenţionale CDMA poate fi atât sincronă cât şi asincronă M Figura 1 8 Alocarea 2W Resurselor în CDMA 1 Resursa de bandă de frecvenţă Resursa de timpT CDMA—sincronă Acest mod este realizabil dacă toate cele M semnale ale 1 utilizatorilor pot fi transmise sincron, de exemplu în legătura descendentă a telefoniei mobile Fiecare utilizator utilizează secvenţa sa unică (numită „semnătură” sau „cheie”) Fiecare semnătură ocupă toată resursa timp-frecvenţă disponibilă Nu există o distribuire a resurselor de timp sau frecvenţă între utilizatori Toate semnăturile sunt sincrone, adică staţia de bază le generează în acord cu un tact propriu Într-o bandă W şi într-un timp T pot opera peste M=W⋅T utilizatori fără interferenţă mutuală, deoarece există M=W⋅T semnale ortogonale, ce joacă rolul de semnături Exemplu: legătura descendentă pentru telefonia mobilă IS-95 dispune de 64 de funcţii Walsh cu rol de semnături pentru 64 de canale ale utilizatorilor Este posibilă creşterea numărului utilizatorilor peste un W⋅T dat, permiţând un anume nivel al interferenţei mutuale Atunci, acest mod devine oarecum asemănător cu CDMA asincron Exemplu: legătura descendentă pentru 3G cdma 2000 CDMA—asincronă se utilizează atunci când este imposibilă sincronizarea tuturor semnăturilor, de exemplu în legătura ascendentă (de la MS către BS) din telefonia mobilă 1 Legătura descendentă, sau inversă (down link) este legătura dinspre o staţie de bază (BS) către o staţie mobilă (MS), sau canalul „unul către mai mulţi” Legătura ascendentă, sau directă (up link) este legătura dinspre o staţie mobilă către o staţie de bază, sau canalul „mai mulţi către unul” 9 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta Diferenţa de timp dintre semnalele utilizatorilor este necontrolabilă şi aleatoare, astfel că semnăturile (spre deosebire de CDMA sincronă) nu pot rămâne ortogonale şi interferenţa mutuală devine iminentă Nivelul ei defineşte numărul de utilizatori Semnăturile trebuie să fie alese în aşa fel încât, dacă un receptor este acordat pe o semnătură, toate celelalte apar ca şi zgomotul pentru acel receptor În acest fel, ţinând cont şi de interferenţa mutuală, RSZ al filtrului adaptat devine 2 ξ 2E2⋅ = = rξ 2 1⋅⋅⋅+ξ T)W2/(1)-(M N+ i0N TW2⋅⋅ (1 12) ≈ M unde Ps / W (1 13) Ni = (M-1) ⋅ 2 = 2⋅E/N0 este este densitatea spectrală de putere a interferenţei totale de la alţi utilizatori iar ξ 2 fixat, numărul maxim de utilizatori este: RSZ în lipsa interferenţei Pentru un ξ TW2⋅⋅ M ≈ (1 14) 2 ξ Comparaţ ie între FDMA, TDMA şi CDMA Pentru FDMA, TDMA şi CDMA sincronă 21-2 Dacă Pb 5,4 (7,3dB), ceea ce conduce la M 1 Numărul de polinoame primitive de grad m este dat prin formula: 20 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta Jm − − i1P12 (2 8) NP(m) = ∏ ⋅ = 1iPm i m {Pi, i = 1, 2, …, J} este descompunerea în numere prime pentru 2 –1, adică J me i 2 – 1 = ∏ iP, =1i unde ei este un întreg În Anexa A este prezentată o listă cu polinoamele primitive, până la gradul 13 Valori tipice de interes sunt pentru m între 10 şi 50 2 2 2 Proprietăţile secvenţelor M În continuare se va demonstra că secvenţele M satisfac îndeaproape proprietăţile aleatoare ale secvenţelor aleatoare binare după cum au fost anunţate în secţiunea 2 1 De remarcat faptul că parametrii procesului de generare [coeficienţii ci ai polinomului generator g(D)] sunt determinişti Singurii parametri aleatori sunt cei m termeni ai vectorului iniţial sau, echivalent, timpul de deplasare R 1: Proprietatea de simetrie Examinăm prima celulă a registrului din Figura 2 1, atunci când întreaga secvenţă este deplasată prin el Vom urmări acest ultim bit (cel mai din dreapta) al vectorului m-dimensional la fiecare ciclu de tact Urmărind conţinutul registrului de deplasare atunci când el generează o secvenţă M, se va vedea că registrul trece prin toate cele m 2-1 stă rile posibile, exceptând-o pe cea nulă Dar dacă am include şi pe cea nulă, am m avea cu siguranţă o secvenţă echilibrată: dintre cei 2 vectori binari, jumă tate sunt pari (au 0 pe ultima poziţie) şi jumătate sunt impari (au 1 pe ultima poziţie) Deoarece am mm inclus starea nulă, echilibrul diferă cu 1 faţă de 2 Aş adar, din 2-1 termeni ai m-1m-1 secvenţei M, 2 sunt unu-uri ş i 2-1 sunt zerouri Dacă presupunem vectorul iniţial aleatoriu, sau echivalent, presupunem un moment de start ales aleatoriu, probabilitatea ca la un anume tact ieşirea registrului să fie zero, respectiv unu, este: 1-m 11 12−  Pm(0) = =  −⋅ 1, m P2 12−  1-m 11 2  Pm(1) = =  +⋅1 m P2 12−  -3-9 Astfel nesimetria este 1/P Pentru m = 10, 30 şi 50, 1/P este aproximativ 10, 10, ş i -15 respectiv 10 R 2: Proprietatea de lungime de fugă Considerăm toate conţinuturile de forma x x …x 0 a1 a2 … ak 0 ş i x x …x 1 b1 b2 … bk 1 21 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta a de unu-uri, „b1 b2 … bk” este secvenţa nulă iar unde „a1 a2 … ak” este secvenţ prefixul „x x …x” poate să fie orice secvenţă binară de lungime m–2–k Deoarece, aşa m cum s-a demonstrat anterior, în decursul unei perioade de 2 –1 tacte, registrul va conţine toate secvenţele posibile de lungime m (exceptând-o pe cea nulă), rezultă un m-2+k număr de 2⋅2 secvenţ e având una dintre formele invocate Acesta reprezintă de fapt numărul de lungimi de fugă de mărime k În plus, numărul total de lungimi de 2-m2-m 1-k-mim-1 fugă de mărime k m R 3: Proprietatea de deplasare şi adunare Considerăm pentru o secvenţă M oricare două perioade de tact deplasate Dacă m luăm întreaga secvenţă de lungime P = 2 – 1 ş i o deplasăm cu un număr arbitrar de perioade τ 20), o secvenţă M este indiscutabil o secvenţă binară Bernoulli sau „coin-flipping”, cel puţin cu respectarea proprietăţilor (R 1) la (R 3), atâta timp cât vectorul iniţial sau timpul sunt alese arbitrar 2 3 Coduri utilizate pentru împrăştierea spectrului Codurile care se pot găsi în sistemele DS practice sunt: coduri Walsh-Hadamard, secvenţe M, coduri Gold şi coduri Kasami Aceste seturi de coduri pot fi în mare împărţite în două clase: coduri ortogonale şi coduri neortogonale Secvenţele Walsh fac parte din prima categorie, în timp ce celălalt grup conţine aşa numitele secvenţe generate cu registre de deplasare 2 3 1 Codurile Walsh Hadamard Secvenţele Walsh Hadamard au avantajul de a fi ortogonale, şi în acest fel se pot elimina interferenţele de multi-acces Există totuşi şi câteva neajunsuri: Codurile nu au un singur maxim al autocorelaţiei • • Împrăştierea nu este pe întreaga lăţime de bandă, energia fiind dispersată pe un număr de frecvenţe discrete (vezi Figura 2 2) 23 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta • Deşi intercorelarea de-a lungul întregii secvenţe este zero, intercorelarea calculată pentru o porţiune din secvenţe este nenulă Consecinţa este că se pierde avantajul de a folosi coduri ortogonale Ortogonalitatea este de asemenea afectată de efectele canalui asupra semnalului ca de • exemplu apariţia căilor multiple În sistemele practice este utilizată egalizarea pentru a recupera semnalul iniţial Aceste neajunsuri fac ca secvenţele Walsh să nu poată fi utilizate în sistemele noncelulare Sistemele în care sunt folosite secvenţele Walsh sunt sistemele celulare CDMA Spre exemplu, sistemul IS-95 foloseşte o combinaţie între o secvenţă Walsh şi o secvenţă M pentru a face posibilă sincronizarea Secvenţele Walsh se obţin prin eşantionarea funcţiilor cu acelaşi nume Funcţiile Walsh formează un set complet ortonormat de funcţii rectangulare Există trei grupuri de funcţii Walsh care diferă între ele numai prin modul de ordonare: – ordonare pe bază de secvenţă (Walsh); – ordonare naturală (Hadamard); – ordonare diadică (Paley) trată de ordin k, ce se poate Secvenţele Walsh sunt linii ale matricii Hadamard, Hk–matrice pă obţine recursiv: ++11   22HH H1 = [+1] H2 = H3 =    −+11H-H  22    1-n1-nHH Hn = (2 11)   1-n1-nH-H   –––secvenţă Walsh - - -secvenţă M Figura 2 2 Comparaţie în domeniul frecvenţă între o secvenţă Walsh şi una M, 24 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta Coduri ortogonale de lungime variabilă Unul dintre obiectivele CDMA este îmbună tăţirea eficienţei comunicaţiilor multimedia Astfel, datorită varietăţii serviciilor, ratele de bit necesare variază de la valori mici la foarte mari Deoarece banda semnalelor de împrăştiere este aceeaşi pentru toţi utilizatorii, diversele rate de transmisie necesită factori de împrăştiere (sau factori expansiune i în canalele fizice Considerăm că fiecare bit al celei mai mici rate de bit de bandă, β e) diferiţ n tiat printr-un cod de lungime N=2 Deoarece durata bitului pentru rata (Rmin) este împrăş Rmin este jumătate din durata bitului ratei minime, este necesar un cod de împrăştiere de 2⋅ n-1kn-k 2⋅Rmin este necesar un cod de lungime 2 lungime N/2=2 În general, pentru o rată Domeniul lungimii codului depinde aşadar de ratele de bit, maximă şi minimă suportate de sistem, şi de lăţimea de bandă O metodă de obţinere a codurilor ortogonale de lungime variabilă, ce păstrează ortogonalitatea între diferitele rate şi factori de împrăştiere, bazată pe transformata Hadamard modificată, este prezentată în rime NxN ce denotă setul de N coduri binare de împrăştiere de Fie CN o matrice de mă n CN(n)}, unde CN(n) este vectorul linie de N=2 elemente; Această N chipuri lungime, { n=1,…,N matrice este generată din CN/2 astfel: )1(C )1(C)1(C   N⋅N/2N/2   )2(C⋅ N)1(C)1(CN/2N/2     )3(C⋅ N)2(C)2(CN/2N/2 C =  =   N (2 12) MM     1)-N(C⋅ NN/2)(CN/2)(CN/2N/2   N)(C N/2)(CN/2)(C ⋅   NN/2N/2   Aşadar, aceste coduri ortogonale de lungime variabilă pot fi generate recursiv utilizând , la structura arborescentă prezentată în Figgura 2 3 Pornind de la C1(1) =1, se generează kk tiere de lungime 2 chipuri nivelul k, un set de 2 coduri de împrăş C8(1) C4(1) = {1,1,1,1} C8(2) C2(1) = {1,1} C8(3) C4(2) = {1,1,-1,-1} C8(4) C1(1) = {1} C8(5) C(3) = {1,-1,1,-1} 4 C8(6) C2(2) = {1,-1} C(7) 8 C4(4) = {1,-1,-1,1} C8(8) e =1 βe =2βe =4βe =8 β Figura 2 3 Arborele utilizat la generarea codurilor ortogonale de lungime variabilă 25 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta Codurile din acelaşi nivel constituie setul de funcţii Walsh şi, ca atare, sunt ortogonale Mai mult, oricare două coduri ale diferitelor niveluri sunt de asemenea ortogonale, exceptând cazul în care unul dintre cele două coduri este părinte al celuilalt cod: de exemplu C16(2), C8(1), C4(1) ş i C2(1) sunt părinţi pentru C32(3) şi, ca atare, nu sunt ortogonali cu acesta Cu alte cuvinte, un cod poate fi utilizat într-un canal dacă şi numai dacă nu este utilizat nici un alt cod al căii codului specificat sau al altei căi ce conduce la aceeaşi rădăcină Aceasta înseamnă că numărul codurilor disponibile nu este fix, ci depinde de ratele de bit şi factorii de împrăştiere ai canalului fizic Această restricţie este impusă cu scopul de a menţine ortogonalitatea 2 3 2 Secvenţe M Secvenţele M, , sunt secvenţe generate prin registru de deplasare (secvenţe MLSR) Secvenţele M nu sunt ortogonale, dar au un maxim îngust al autocorelaţiei Aşa cum s-a arătat deja, aceste coduri pot fi generate utilizând un registru de deplasare Obţinerea de secvenţe de secvenţe M este condiţionată de utilizarea conexiunilor buclei de reacţie în acord cu coeficienţii unui polinom primitiv de grad m În Anexa A sunt prezentate polinoamele m primitive până la gradul m=13 Lungimea secvenţei M generate este egală cu 2-1 Numă rul de coduri posibile este dat de numărul de polinoame primitive existente pentru respectivul m, de asemenea prezentat în Anexa A Aşa cum s-a arătat în paragraful 2 2 2, secvenţele M au câteva proprietăţi speciale: • Secvenţele M sunt cvasi-echilibrate: numărul de 1 este cu 1 mai mare decât numărul de 0 2 • Spectrul unei secvenţe M are o anvelopă de forma „sinc” În Figura 2 2 este prezentat, , spectrul unei secvenţe Walsh de lungime 64 şi al unei secvenţe M de lungime 63 Ambele secvenţe sunt (aproape) de aceeaşi putere Figura arată că aplicarea unei secvenţe M distribuie mai bine puterea pe întregul domeniu de frecvenţe decât secvenţa Walsh Proprietatea „deplasează-şi-adună” poate fi formulată după cum urmează: • kij T u = T u + T u (2 13a) Aici u este o secvenţă M Combinând două deplasări ale acestei secvenţe (deplasări relative i şi j) obţinem secvenţa M, cu o altă deplasare relativă • Funcţia de autocorelare este bivalentă: Nτ Nk ⋅=  () Rτ =  u (2 13b) 1τ Nk ⋅≠−  unde k este o valoare întreagă, iar τ este deplasarea relativă • Nu există o formulă generală pentru intercorelarea a două secvenţe M, putând fi formulate doar câteva reguli • O aşa numită „pereche preferată” este o combinaţie de secvenţe M pentru care [(m+1)/2][(m+1)/2] intercorelarea duce la doar 3 valori diferite: -1, -2 ş i 2 –2 Nu există perechi preferate pentru secvenţele M cu o lungime de 4⋅k unde k este un întreg 2 3 3 Secvenţe Gold şi Kasami În unele aplicaţii proprietăţile intercorelaţiei secvenţelor PN sunt la fel de importante ca şi proprietăţile de autocorelaţie Spre exemplu, în CDMA fiecărui utilizator îi este atribuită 26 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta o secvenţă PN particulară Ideal, secvenţele PN dintre utilizatori ar fi mutual ortogonale însă, secvenţele PN utilizate în practică prezintă o oarecare corelaţie Funcţia de intercorelaţie periodică dintre orice perechi de secvenţe M de aceeaşi perioadă poate avea maximuri mari Tabelul 2 4 prezintă valorile maximurilor Φ max pentru intercorelaţia dintre perechi de secvenţe M pentru 3≤m≤12 Este prezentat, de asemenea, în m Tabelul 2 4 numărul de secvenţe M de lungime n = 2-1 pentru 3≤ m≤12 Aşa cum se poate vedea din acest tabel, numărul de secvenţe de lungime n creşte rapid cu m De asemenea observăm că, pentru multe secvenţe, mărimea maximului Φ iei este un max a intercorelaţ procent mare din valoarea funcţiei de autocorelaţie, Φ 0 Astfel de valori mari ale intercorelaţiei sunt inacceptabile în CDMA Chiar dacă este posibil să se selecteze un mic subset de secvenţe M care au valori relativ mici ale maximului intercorelaţiei, numărul de secvenţe al subsetului este uzul prea mic pentru aplicaţii CDMA Secvenţe PN cu proprietăţi ale intercorelaţiei periodice superioare secvenţelor M au fost date de Gold (1967, 1968) şi Kasami (1966) Ele sunt derivate din secvenţele M aşa cum se va descrie în cele ce urmează Gold şi Kasami au dovedit că există perechi de secvenţe M, de lungime n, ce dau trei valori pentru funcţia de intercorelaţie {-1, -t(m), t(m)-2} unde 2/1m () +  12 m impar =+ t(m) =  (2 14) () +2/2m 12 m par =+  6 De exemplu, dacă m = 10, t(10) = 2 + 1 = 65 iar cele trei posibile valori pentru funcţ ia de intercorelaţie sunt {-1, -65, 63} Astfel maximul intercorelaţiei pentru perechile de secvenţe M este 65 în vreme ce maximul familiei de 60 de secvenţe posibile generate printr- un registru de deplasare cu 10 stări, cu diferite conexiuni de reacţie, este Φ max = 383, de aproape 6 ori diferenţa în valoarea maximă i p2(x), polinoame primitive de grad m ce Două secvenţe M, generate de p1(x) ş t(m)m constituie polinoamele minimale pentru elementele α şi α ale câmpului Galois GF(2), a căror funcţie de intercorelaţie ia trei valori posibile {-1, -t(m), t(m)-2} se numesc secvenţ e ele preferate posibile pentru m ≤12 preferate În Anexa B sunt prezentate secvenţ m m n=2 -1 Numărul de Maximul Φ max/Φ0 t(m) T(m)/ Φ0 secvenţe M intercorelaţiei Φ max 3 7 2 5 0,71 5 0,71 4 15 2 9 0,6 9 0,6 5 31 6 11 0,35 9 0,29 6 63 6 23 0,36 17 0,27 7 127 18 41 0,32 17 0,13 8 255 16 95 0,37 33 0,13 9 511 48 113 0,22 33 0,06 10 1023 60 383 0,37 65 0,06 11 2047 176 287 0,14 65 0,03 12 4095 144 1407 0,34 129 0,03 Tabelul 2 4 Vârful intercorelaţiei secvenţelor M şi Gold 27 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta i b = [b1 b2 … bn], Pentru o pereche de secvenţe preferate, să spunem a = [a1 a2 … an] ş construim un set de secvenţe de lungime n prin sumarea modulo 2 a secvenţei a cu versiunea 1 permutată ciclic a secvenţei b, sau viceversa Atunci obţ inem n noi secvenţe periodice de m perioadă n=2-1 Putem include, de asemenea, secvenţ ele originale a şi b şi, astfel, avem în total n+2 secvenţe Cele n+2 secvenţe construite în această manieră se numesc secvenţ e Gold Cu excepţia secvenţelor a ş i b, secvenţele Gold nu sunt secvenţe M (nu au perioadă maximă) Prin urmare funcţiile lor de autocorelaţie nu au doar 2 valori Gold (1968) a arătat că funcţiile de autocorelaţie pentru orice pereche din setul de n + 2 secvenţe Gold are trei valori posibile {-1, -t(m), t(m)-2}, unde t(m) este dat prin (2 14) Similar, maximul secund al funcţiei de autocorelaţie pentru o secvenţă Gold ia valori din setul {-1, -t(m), t(m)-2} Rezultă că valorile vârfului (secund) al funcţiei de autocorelaţie este mărginit superior de t(m) Valorile maximului secund al autocorelaţiei şi maximului intercorelaţiei, adică t(m), pentru secvenţele Gold sunt listate în Tabelul 2 3 Sunt de asemenea listate valorile normalizate pentru Φ(0) Este important de comparat valoarea maximului intercorelaţiei pentru secvenţele Gold cu marginea inferioară a intercorelaţiei dintre orice pereche de secvenţe binare de perioadă n Marginea inferioară dedusă de Welch (1974) pentru Φ max este: 1-M n⋅ (2 15) Φ max≥ 1-nM⋅ ele Gold, care, pentru valori mari ale lui n, este bine aproximată ca fiind n Pentru secvenţ mm/2 n= 2 –1 ş i, de aceea, marginea inferioară este Φ ≈ 2 Marginea inferioară Welch este max e Gold i cu 2 pentru m par decât Φ mai mică cu 2pentru m impar ş max = t(m), pentru secvenţ O procedură similară cu cea folosită pentru generarea secvenţelor Gold va genera setul m/2m e binare de perioadă n = 2 –1, unde m este par În această procedură, de N = 2 secvenţ m/2 pornim cu o secvenţă M, a, ş i formăm o secvenţă binară b luând tot al 2+1 –lea bit din a m/2 Adică secvenţa b este formată decimând a prin 2+1 Se poate verifica că rezultatul, b, este m/2 m=10, perioada lui a este n = 1023 o secvenţă periodică cu perioada 2-1 De exemplu, dacă iar perioada lui b este 31 Astfel, dacă observăm 1023 biţi pentru secvenţa b, vom vedea 33 de m i pentru secvenţele a şi b, formăm repetiţii de 31 de secvenţe de biţi Luând acum n = 2-1 biţ m/2 un nou set de secvenţe prin adunarea, modulo 2, a lui a cu b ş i tuturor celor 2-2 secvenţe m binare de lungime n = 2-1 Acestea sunt denumite secvenţ e Kasami Funcţiile de m/2 autocorelaţie şi intercorelaţie pentru aceste secvenţe iau valori din tripletul {-1, -(2+1), m/2 2-1 } Astfel, valorile maximului intercorelaţiei pentru orice pereche de secvenţe ale setului este: m/2 Φ max = 2 + 1 (2 16) m/2 Aceste valori pentru Φ Welch pentru un set de 2 secvenţe max satisfac limitarea inferioară m de lungime n = 2 –1 Din această cauză secvenţele Kasami sunt optimale 2 3 4 Concluzii Soluţii alternative Alegerea codurilor (secvenţ elor) PN depinde de aplicaţie Fiecare dintre codurile PN prezentate oferă anumite avantaje Secvenţele Walsh-Hadamard fiind ortogonale, vor fi utile 1 O metodă de generare a celor n secvenţe este printr-un registru de deplasare de lungime 2⋅m cu conexiunile reacţiei specificate prin polinomul h(p) = g1(p)⋅ g2(p), unde g1(p) şi g2(p) sunt polinoamele care specifică conexiunile reacţiei registrelor de deplasare cu m stări ce generează secvenţele M, a şi b 28 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta în sistemele de transmisie cu spectru împrăştiat cu CDMA unde se poate realiza sincronizarea tuturor convorbirilor Secvenţele M sunt uşor de generat, însă nu au proprietăţi bune de intercorelare Nu vor putea fi folosite în sisteme multiacces Astfel de proprietăţi de intercorelaţie au secvenţele Gold sau Kasami Însă, în multe sisteme CDMA practice durata bitului nu acoperă decât o fracţie din secvenţa periodică Într-un astfel de caz, este importantă intercorelaţia periodică parţială dintre două secvenţe În este prezentat modul de construcţie a unui set de secvenţe Gold optimal din punct de vedere al intercorelaţiei parţiale Secvenţe PN pot fi privite ca şi cuvinte de cod protector, ceea ce le face utile la recepţie Astfel, codurile Hadamard, aşa cum se va prezenta în Cap 4, sunt cuvinte de cod Reed-Muller, iar secvenţele M sunt de fapt cuvinte de cod ciclic, fiind generate acelaşi procedeu (Cap 5) Acest fapt permite utilizarea tehnicilor de decodare ale codurilor protectoare şi pentru secvenţe PN, Selecţia secvenţelor preferate devine o problemă dificilă pentru m mare În este prezentată o metodă de aflare a unor astfel de perechi de secvenţe Spre exemplu două secvenţe preferate pentru m =127 sunt: 127 x + x +1, şi 12711399988885848071695957 x + x + x + x + x + x + x + x + x + x + x + x + 555144434241403735322928 + x + x + x + x + x + x + x + x + x + x + x + x + 26252423201911105 + x + x + x + x + x + x + x + x + x + x + 1 În sistemele cu spectru împrăştiat sunt utilizate şi secvenţe neperiodice, numite haotice Ele sunt generate de sisteme neliniare Aceste secvenţe haotice oferă, în comparaţie cu primele, avantajul unor mai bune proprietăţi aleatoare (de zgomot) şi posibilitatea unei sincronizări mai precise la recepţie, O altă soluţie alternativă secvenţelor „clasice”, o reprezintă secvenţele nebinare Acestea posedă proprietăţi superioare din punct de vedere al intercorelaţiei, 29 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta 3 Tehnici adiţionale de îmbunătăţire a performanţelor sistemelor cu spectru împrăştiat Prin natura sa, tehnica împrăştierii spectrului oferă o resursă timp-frecvenţă considerabil mai mare decât cea necesară transmisiei informaţiei pentru oricare utilizator în parte Aceasta se reflectă prin câştigul de procesare, Gp, mare, sau prin raportul W/R (lăţ ime de bandă–rată de transmisie) mare Acest exces de redundanţă poate fi exploatat pentru a îmbunătăţi performanţa, fără a compromite alte avantaje ale câştigului mare de procesare În acest capitol sunt prezentate două tehnici de procesare ce aduc îmbunătăţiri sistemului de transmisie: întreţeserea (interleavingul) şi codarea pentru corecţia erorilor 3 1 Întreţeserea Avantajul imediat al excesului de redundanţă este că determină independenţa ieşirilor canalului Cu cât e mai mare numărul de componente de cale independente disponibile în prezenţa fadingului, cu atât e mai bună performanţa (vezi paragraful 1 2 punctul g) Dispunem de L căi şi de N chipuri per simbol Problema e că chipurile succesive nu pot fi privite ca independente Este posibil, însă, să reordonăm chipurile astfel încât cele N chipuri ce provin de la acelaşi simbol să nu mai fie transmise succesiv Mai exact, ele sunt transmise la intervale suficient de largi astfel că fading-ul va conduce la amplitudine şi fază Ieşire către canal Intrare dinspre canal J coIeş Intrare loane J coloaneire x1 x2 x3 x4 xJ x1 x2 x3 x4 xJ x1 x2 x3 xJ+ xJ+ 11 x2J x2J+1 +1 I rânduri I rânduri x(I-1)J+1 xIJx(I-1)J+1 xIJ Interleaver Deinterleaver IE pe COLOANE ÎNSCRIERE pe RÂNDURIRECEPŢ CITIRE pe COLOANE CITIRE pe RÂNDURI Secvenţa întreţesută: x1, xJ+1, x2J+1, … ,x(I-1)J+1, x2, xJ+2 Oricare două simboluri aflate iniţial la mai puţin de J poziţii unul de celălalt vor fi depărtate la cel puţin I poziţii Figura 3 1 Întreţesere/deîntreţesere bloc 30 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta 1 independente pentru fiecare dintre cele N chip-uri Procedeul de reordonare necesar pentru a dobândi diversitate temporală este numit întreţesere (interleaving) şi poate fi realizat în mai multe feluri Două procedee sunt întreţeserea „bloc” şi întreţeserea „convoluţională” (Figura 3 1 şi 3 2) După recepţie, chipurile demodulate de la ieşire sunt din nou reordonate pentru a le pune la loc în ordinea originală prin procedeul invers numit deîntreţesere (deinterleaving) Acest proces evident introduce întârziere între generarea datelor digitale şi livrarea lor la utilizatorul receptor Întârzierea produsă prin întreţeserea bloc este de aproximativ 2⋅I⋅J⋅T (T este durata unui simbol) Aşa cum se poate vedea din cele două figuri, întârzierea produsă de întreţeserea convoluţională este aproximativ jumătate din întârzierea produsă de întreţeserea bloc: J⋅(I–1)⋅T, Ordinea la intrarea interliverului şi la ieşirea deinterliverului: … xi, xi+1, xi+2, xi+3, … Ordinea la ieşirea interliverului şi la intrarea deinterliverului: I xi, xi-J, xi-2J, …, xi-(I-1)J, xi+I, xi+I-J, …, xi+1 J Figura 3 2 Întreţesere/deîntreţesere convoluţională Aceasta este doar un alt exemplu al faptului cunoscut că, în canalele cu fading, diversitatea poate îmbunătăţi considerabil performanţa Această diversitate poate fi realizată prin separare spaţială, cu antene multiple sau prin căi multiple produse natural, sau prin separare temporală cu ajutorul procesului de întreţesere descris anterior Preţul ultimei metode este întârzierea Cu cât procesul de fading este mai lent, cu atât este necesară o întreţesere pe o întindere mai mare şi de aceea o întârziere mai mare pentru realizarea independenţei prin diversitate temporală 1 Întreţeserea poate fi (şi uzual este) făcută nu pe chipuri individuale, ci pe sub-secvenţe de N chipuri, înainte de împrăştierea pseudo-aleatoare Aceasta reduce necesarul de memorie, deoarece toate chipurile din fiecare sub- secvenţă, care au fost ţinute apropiate după întreţesere, au acelaşi semn Mai mult, codând, sub-secvenţa constă în general dintr-un cuvânt de cod Întreţeserea simbolurilor codate conferă mai mare avantaj fără memorie excesivă şi necesităţi de procesare pentru întreţeserea de chip 31 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta Întreţeserea este utilizată şi în combinaţie cu concatenarea codurilor, aşa cum se va prezenta în paragraful 3 3 3 Performanţele codurilor respective depind în mod evident şi de dispozitivul de întreţesere utilizat Din nefericire, performanţele întreţeserii depind de caracteristicile canalului încât nu există soluţii universal valabile în ceea ce priveşte dispozitivul de întreţesere, 3 2 Codarea corectoare de erori Metoda efectiv universală de exploatare a redundanţei este codarea corectoare de erori (FEC–forward error-correcting) FEC îmbunătăţeşte performanţa pentru canalele cu amplitudine şi fază fixe la fel de bine ca şi pentru canalele cu fading În capitolele următoare vor fi prezentate câteva dintre cele mai utilizare coduri în sistemele cu spectru împrăştiat 3 2 1 Decodare soft / hard Indiferent de codul utilizat, structura unui receptor conţine un bloc de decizie Acest bloc reprezintă un comparator ce compară doi termeni pentru a hotărâ valoarea fiecărui bit de informaţie din secvenţa recepţionată: D1 Λ ; (3 1) i µ D0 celui de-al Dacă primul termen, Λ i, numit raport de plauzibilitate, este de valoare superioară doilea, µ, numit pragul deciziei, atunci valoarea bitului în cauză este “1” În mod curent Λ i depinde de semnalul recepţionat şi este un raport construit pentru fiecare simbol binar în parte, iar µ este constant pentru o aplicaţie dată Dacă decizia se face înainte de decodare, decodarea se numeşte hard Dacă, vice versa, decizia se execută după decodare, decodarea se numeşte soft În acest din urmă caz, blocul de decizie este parte din cecodor (se află la finele acestuia) În cazul decodării hard intrarea decodorului este o secvenţă binară Decodorul o analizează şi, eventual, constatând prezenţa erorilor, corectează inversând valoarea binară a biţilor presupuşi eronaţi În cazul decodării soft intrarea decodorului este un semnal numeric multinivel, rezultat prin eşantionarea (eventual şi cuantizarea valorilor eşantioanelor) semnalului recepţionat, şi/sau, funcţie de aplicaţie, prin demodulare Analizând semnalul numeric de la intrarea sa, decodorul construieşte pentru fiecare bit termenul Λ i-l compară cu pragul deciziei, µ i ş (execută decizia) Utilizarea decodării soft creşte câştigul de codare Decodarea soft poate fi utilizată în combinaţie cu concatenarea codurilor 3 2 2 Câştigul de codare Câştigul de codare, pentru o aplicaţie dată, se defineşte ca şi raportul între puterea semnalui, Ps, necesară pentru a se obţine o anumită rată, impusă, a erorii, pentru transmisia necodată şi puterea semnalului, Psc, necesară pentru a obţine aceeaşi rată a erorii în cazul transmisiei codate: Gc = 10⋅ lg (Ps/Psc) [dB] (3 2) 32 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta Presupunând acelaşi nivel de zgomot (în cazurile “codat” şi “necodat”), rezultă că Gc se poate defini şi ca diferenţa între rapoartele semnal per zgomot necesare obţinerii ratei de eroare impuse, în cele două cazuri, necodat şi codat:  0sNPsPscP  Gc = 10⋅ lg= 10⋅lg –10⋅lg= ξ – ξ (3 3) ⋅ ssc  sc0PN0N0N  3 2 3 Criteriul MAP Conform criteriul MAP (Maximum Aposteriori), prin procesul de decodare se caută acea secvenţă emisibilă vk, ce maximizează probabilitatea condiţionată aposteriori p(vj/r), peste toate secvenţele emisibile {vj} j=1÷N Indiferent de metoda de decodare propriu-zisă, algoritmul decodării urmăreşte, principial, aflarea acelei secvenţe de cod (cuvânt de cod), vk, din mulţ imea secvenţelor posibil a fi emise, V = {vj} probabilitatea aposteriori: j=1÷N, ce maximizează P(k) = p(vk/r) (3 4) unde r este secvenţa recepţionată curentă Cu alte cuvinte, criteriul de selecţie MAP spune astfel: „dacă s-a recepţionat r, alege drept cuvânt emis pe vk, unde p(vk/r) > p(vj/r) ∀ vj = cuvânt emisibil (∈V), j ≠ k” (3 5) Dacă admitem ipoteza independenţei transmiterii biţilor (canal fără memorie), şi în plus po 1,trebuie să fie peste valoarea q Spre exemplu, dorind a construi un cod corector de t = 2 erori, cantitatea de informaţie necesară a fi furnizată de biţii de control, în număr de m, trebuie să fie: n/2] ≈ 2⋅q –1, (5 51) m ≥ ip1 + ip2 = log2 [1 + n + (n +1)⋅ adică o informaţie aproape dublă (şi implicit un număr dublu de biţi de control) faţă de cazul corecţiei unei erori q o Deoarece precizarea un caracter din cele n+1 ale câmpului Galois GF(2) reprezintă informaţie de q biţi, rezultă posibilitatea de a crea „legături” între biţii de informaţie şi cei de control pentru un cuvânt de cod BCH, construind polinomul generator, g(x), astfel încât să aibă rădăcini t elemente ale câmpului respectiv În felul acesta, deoarece orice cuvânt de cod este multiplul lui g(x), cele t rădăcini ale lui g sunt rădăcini şi pentru orice v Această proprietate a cuvântului de cod furnizează informaţia necesară şi suficientă pentru ca decodorul să afle poziţia a t erori Obs Cele t rădăcini invocate în construcţia polinomului g(x) nu pot fi alese oricare dintre cele n+1 Asta deoarece g(x) trebuie să fie un polinom cu coeficienţi binari g(x) este de fapt un produs de polinoame minimale, a căror rădăcini sunt elementele respective Polinomul minimal mk(x) ce are ră dăcină pe tk, este prin definiţie polinomul cu coeficienţi binari, de gradul cel mai mic posibil, ce are ca rădăcină pe tk (pentru mai multe detalii vezi ) 59 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta 5 4 2 Decodarea codurilor BCH Aşa cum s-a arătat în paragraful precedent decodarea codurilor BCH nu este posibilă prin acelaşi procedeu ca şi la codurile ciclice corectoare de o eroare Decodarea codurilor BCH este prezentată algoritmic, în cele ce urmează: Pasul 1 Se verifică dacă există erori în cuvântul recepţionat În acest scop se calculează coeficienţii sindrom: 2j-1 ) cu j = 1÷t (5 52) S2j-1 = w(α 2j-1 , cu j = 1÷t, sunt cele t unde w(x) este polinomul ataşat cuvântului recepţionat, iar α rădăcini ale lui g(x) că s-a recepţionat un cuvânt corect şi Dacă S2j-1 = 0 pentru orice j, se concluzionează se trece la Pasul 5 Dacă există S2j-1 ≠ 0, atunci: Pasul 2 Se calculează coeficienţii polinomului: σ(x) = (x + X1)⋅ (x + X2) … (x + Xt) = tt-1 = x + σ x + … + σ⋅x +σ (5 53) 1⋅t-1t ji sunt locatorii erorilor, ij precizând poziţia erorii Coeficienţii σ la σ se unde Xj = α 1t calculează funcţie de coeficienţii sindrom după cum urmează: 2j-1 –ţinând cont de faptul că w = v + ε şi că v(α ) = 0, ∀j = 1÷t, rezultă 2j-12j-1 S2j-1 = w(α ) = ε(α) = tt 1-j2 i1-2j j = () = ∑t (5 54) αX cu j = 1÷ ∑ i =1i=1i –pe de altă parte, înlocuind x cu Xj în (5 53) se poate forma sistemul: t i-tt ⋅Xσ= X cu j = 1÷t (5 55) ∑ jjj = 1i p –înmulţind pe rând ecuaţiile cu X, p = 1÷ i sumându-le obţinem sistemul jt, ş  σSt + σSt-1 + … + σS1 = St+1 1⋅2⋅t⋅  σ⋅St+1 + σ⋅St + … + σ⋅S2 = St+2 (5 56) 12t ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅   σ⋅S2t-1 + σ⋅S2t-2 + … + σ⋅St = S2t 12t Dintre cei 2⋅ t coeficienţi sindrom, Sj, cei pari se află cu relaţia 2 S2t = (St) (5 57) care este o consecinţă a faptului că se operează în câmpul binar 60 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta Pasul 3 Cunoscând coeficienţii σ, cu j = 1÷t, se pot afla locatorii erorilor, Xj din ecuaţia: j t i- ∑ ⋅Xσ jj= 1 (5 58) 1i = căutând soluţii de forma: ji (5 59) Xj = α n⋅j Ecuaţia (5 58) se poate modifica, ţinând cont de identitatea α ≡1: tttt ki⋅kiin⋅ki)n(⋅ji ⋅−−⋅− ∑ ⋅ασ⋅⋅αασ⋅ασ⋅ασ? (5 60) j= ∑j= ∑j= ∑j=1 1i=1i=1i=1i = Pasul 4 Se inversează biţii corespunzători poziţiilor j pentru care este satisfăcută relaţia (5 60) Pasul 5 Se selectează biţii de informaţie, primii k biţi din cuvântul recepţionat, verificat şi eventual corectat 5 4 3 Codul Golay Codul Golay împreună cu codul Hamming corector de o eroare sunt singurele coduri corectoare de erori perfecte Codul Golay are parametrii: n = 23, k = 12 şi m = 11 Este un cod corector de 3 erori Este un cod perfect fiindcă numărul cuvintelor eroare corectabile este egal cu cel al corectorilor: 11321 12CCC− ++ (5 61) 232323= Polinomul generator poate fi ales dintre: 11106542 g1(x) = x + x + x + x + x + x + 1 119765 g1(x) = x + x + x + x + x + x + 1 Decodarea se poate face pe baza corespondenţei dintre corectori şi cuvinte eroare corectabile, prezentată în Anexa C Spre exemplu fie cuvântul recepţionat: w = 313213058 Decodarea sa presupune: i polinom generator: 1 –calculul restului r, presupunând g1 ca ş 1 1 0 0 1 0 1 1 0 1 0 0 0 1 0 1 1 0 0 0 1 0 1 1 1 0 0 0 1 1 1 0 1 0 1 1 1 0 0 0 1 1 1 0 1 0 1 1 0 0 0 1 0 0 0 0 1 0 0 1 1 0 0 0 0 0 1 0 1 0 1 1 1 0 0 0 1 1 1 0 1 0 1 1 1 0 0 0 0 0 1 0 0 0 1 1 1 0 0 0 1 1 1 0 1 0 1 0 0 1 1 0 0 1 0 0 0 1 = 0 6 2 1 = r 61 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta 2 –găsirea cuvântului eroare corespunzător În acest scop se caută în tabel linia corespun- zătoare valorii „062” şi în această linie termenul al doilea (ce corespunde cifrei „1”) Rezultatul căutării fiind „160602”, rezultă că w are erori pe poziţiile 16, 6 şi 2 a fi 3 –se corectează poziţiile găsite eronate, rezultând v = 3 1 2 2 1 3 4 78, care se verifică cuvânt de cod 5 5 Codurile Reed-Solomon Codurile Reed-Solomon (RS) fac parte din categoria codurilor ciclice, însă sunt coduri nebinare , Spre deosebire de celelalte coduri ciclice, alfabetul codului RS nu este câmpul binar {0, 1}, ci un câmp finit de ordin superior, numit câmp Galois şi care va fi descris în paragraful următor În acest fel, cuvintele codului RS nu sunt secvenţe (succesiuni) de biţi, ci de caractere Aceste caractere pot fi reprezentate, la rândul lor, prin secvenţe binare, însă sunt indivizibile din punct de vedere al codării şi decodării Reed-Solomon Structural, cuvintele de cod RS au aceeaşi alcătuire ca şi cele de cod ciclic: q GF(2, p(x)) j = 0÷n-1 (5 62) u = un-1 un-2 u1 u0 uj ∈ unde: u – cuvântul de cod, format din n caractere; ie, în număr de m; un-1 un-2 uk – caracterele de informaţ r de k; uk-1 uk-2 u0 – caracterele de control, în numă q – ordinul câmpului; p(x) – polinomul generator al câmpului GF Relaţia de codare are aceeaşi formă ca şi la codurile ciclice: kk x/g(x)) (5 63) u(x) = i(x)⋅x + rest (i(x)⋅ unde g(x) este polinomul generator al codului, al cărui construcţie este prezentată în paragraful următor, iar m-1m-2 x + un-2⋅x + + uk+1⋅x + uk (5 64) i(x) = un-1⋅ este polinomul de informaţie Prin relaţia de codare (5 63) se obţine polinomul ataşat cuvântului de cod, polinom ai cărui coeficienţi sunt tocmai caracterele ce alcătuiesc cuvântul de cod dat de (5 62) Relaţia (5 63) indică, deasemenea, că u(x) este un multiplu al lui g(x) Codul RS, având parametrii n, k şi m, construit după relaţia (5 63), este capabil să corecteze un număr ec de caractere eronate, unde: 2⋅ec = k= n-m (5 65) La decodare, spre deosebire de codurile ciclice, într-un cuvânt de cod RS recepţionat, în vederea corecţiei, este necesară atât localizarea erorii, cât şi stabilirea valorii ei q 5 5 1 Câmpul Galois GF(2) i reprezintă Câmpul Galois este generat de un polinom primitiv p(x) de grad q ş mulţimea claselor de resturi rezultate prin împărţirea polinoamelor cu coeficienţi binari la 33 p(x) În figura 5 11 este prezentat câmpul Galois GF(2, p(x)=x+x+1) 62 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta j z α M r(α ) 0 0 000 0 1 1 001 1 2 α 010 α 22 3 α100 α 3 4 α011 α+ 1 42 5 α110 α +α 52 6 α111 α +α+1 62 7 α101 α +1 3 Figura 5 11 Câmpul Galois GF(2) În tabel se disting patru reprezentări distincte ale elementelor câmpului: - zecimal sau octal (coloana “z”); j ”); - monom (coloana “α - matricial (coloana “M”); )”) - polinom-clase de resturi modulo p(x) (coloana “r(α Pe mulţimea claselor de resturi p(x) sunt definite operaţiile de adunare (sumă modulo doi) şi înmulţire Vom exemplifica modul de operare al acestor doi operatori Fie două elemente (clase rest): x a = aq-1 aq-2 a1 a0 = α y (5 66) b = bq-1 bq-2 b1 b0 = α Elementul sumă este: z c = a+b = aq-1+bq-1 aq-2+bq-2 a1+b1 a0+b0 = c (5 67) iar elementul produs: x+y (5 68) d = a⋅b = dq-1 dq-2 d1 d0 = α Trebuie menţionat că: 7n0 = α =α =1 (5 69) α Exemplu: fie cele două elemente a şi b din câmpul Galois prezentat în Figura 5 11: 2 a = 3 = 100 = α 4 b = 5 = 110 = α şi Rezultă că c = 010 = 2 = α 6 d = α = 101 = 7 Pentru a putea urmări înmulţirea utilizând doar reprezentarea zecimală, fără a face apel j , se procedează astfel: la coloana α a +z b ≤ 8 d = a⋅b = a +z b –z 1 dacă a +z b > 8 (5 70) = a +z b –z 8 dacă unde: q “⋅” - operatorul înmulţire în GF(2); 63 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta i “–z” - operaţiile de sumare şi scădere din N-mulţimea numerelor naturale “+z” ş q ii de sumare în În concluzie, operaţia de înmulţire în GF(2) corespunde unei operaţ zecimal De asemenea, pentru implementarea fizică, poate fi înţeleasă ca o operaţie de i “b” au reprezentare numărare de la elementul “a” până la elementul “a +z b –z 1”, unde “a” ş zecimală (prima coloană în Figura 5 11) În contraparte, operaţia de divizare a lui “a” la “b” i totodată unei numărări în va corespunde unei diferenţe “a –z (b –z 1)” sau “a –z (b –z 8)” ş la “a –z (b –z 1)” sau “a –z (b –z 8)” sens descrescător de la “az” până 5 5 2 Polinomul generator, g(x), al codului ia fiecăreia Pentru a se corecta t erori dintr-un cuvânt este necesar a se preciza poziţ precum şi valoarea ei Dacă ne referim la un cuvânt de lungime n, unde: q n = 2 – 1 (5 71) atunci informaţia necesară pentru a preciza un caracter eronat între cele n este conţinută de un q ia (5 50) De asemenea şi valoarea erorii, ε, poate să fie orice caracter din GF(2), relaţ q caracter din GF(2): w = v + ε (5 72) q unde: w –caracterul recepţionat ∈GF(2); q v –caracterul emis ∈GF(2); q ε -valoarea erorii ∈GF(2)\{0} Incluzând şi cazul “eronare a caracterului fictiv”, rezultă că ε poate lua orice valoare din qq 2 valori posibile Informaţia necesară pentru a preciza valoarea ei este identică GF(2), adică cu cea dată de (5 50) În concluzie, pentru fiecare eroare ce se doreşte a fi corectată este necesară o q informaţie egală cu 2q biţi, adică două caractere din GF(2) La t erori sunt necesare 2t caractere (cantitate de informaţie) Obs În fapt condiţia anterioară este una suficientă, cea necesară implică mai puţină informaţie deoarece nu se poate erona un caracter de două ori Cele 2t caractere de informaţie necesare soluţionării problemei corecţiei se află din 2t ecuaţii, care înseamnă tot atâtea legături (proprietăţi) pentru cuvântul recepţionat Aceste 2t proprietăţi pentru cuvintele de cod RS sunt generate prin relaţia de codare (5 63) Prin această relaţie u(x) devine multiplul lui g(x), ceea ce înseamnă că rădăcinile lui g vor fi şi rădăcini pentru u Rezultă necesitatea ca g să aibă 2t rădăcini Aceste rădăcini pot fi oricare dintre cele qq232t dăcinile α, α, α, α, 2 elemente ale câmpului GF(2) Spre exemplu, se pot alege ră datorită simplităţii şi simetriei: 22t ) (x + α) g(x) = (x + α)(x + α kk-1 = x + gk-1x + +g1x + g0 (5 73) Aşadar cuvântul de cod RS, u, rezultat prin codarea cu ajutorul relaţiei (5 63), în care g 232t , α, α sunt rădăcini pentru este dat de (5 73), are proprietatea că elementele α, α polinomul ataşat, u(x) 64 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta 5 5 3 Decodarea codurilor Reed-Solomon Decodarea codurilor RS este în mare parte asemănătoare decodării codurilor BCH, deosebirea constând în necesitatea de a afla pentru fiecare eroare pe lângă poziţie şi valoarea ei De remarcat că un cuvânt de cod RS este format din n caractere, adică q⋅n biţi Astfel o eroare de caracter poate însemna până la q erori de bit În cele ce urmează este descrisă decodarea codurilor RS prin paşi algoritmici: Pasul I Conform relaţiei (5 73) g(x) este un polinom de grad k=2t > 2 Prin construcţie, cuvintele de cod RS au proprietatea că polinoamele ataşate lor sunt divizibile cu g(x), adică 222t ,α,…,α sunt rădăcini atât pentru g(x) cât şi pentru orice elementele câmpului GF (2) α cuvânt de cod u(x): j ) = 0 g(α j ) = 0 j = 1÷2t (5 74) u(α Aceste proprietăţi constituie şi punctul de plecare în decodare Presupunând că V este un cuvânt recepţionat: v=u+ε sau v(x) = u(x) + ε(x) (5 75) vom calcula 2t coeficienţi, numiţi coeficienţi sindrom, Sj, în forma: jjjj ) = u(α) + ε(α) = ε(α) j=1÷2t (5 76) Sj = v(α Dacă nu există erori Sj = 0 În acest caz se trece la pasul VI Evident concluzia poate fi eronată Un exemplu în argumentarea acestei afirmaţii este situaţia: ε = cuvânt de cod Dar în acest caz numărul erorilor depăşeşte puterea de corecţie de t erori Pasul II Dacă există erori în limitele corectabile (numărul erorilor este mai mic sau egal cu t) atunci i sindrom diferiţi de zero Fie cuvântul eroare în forma: există coeficienţ t kr ii (5 77) () ⋅=xxαε ∑ i1 = unde: r (5 78) iYα i= {0, 1, 2, … n-1}, iar: reprezintă valoarea erorii ri ∈ k (5 79) iXα i= {0, 1, 2, …, n-1} Cu aceste notaţii coeficienţii sindrom au reprezintă locatorul erorii ki ∈ expresiile: tt jkji 2t (5 80) () ⋅=⋅=XYYSα ∑∑ iiij , j = 1÷ ii11 == 65 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta Ecuaţiile (5 80) reprezintă un sistem de 2t ecuaţii cu 2t necunoscute: t locatori ai i t valori pentru respectivele erori Yi Rezolvarea acestui sistem de ecuaţii se va erorilor Xi ş face în mai multe etape La pasul prezent se vor calcula coeficienţii polinomului σ(x) ai cărui rădăcini sunt locatorii erorilor: t − tt1 ( ) (5 81) () +⋅++⋅+=+=xxxXxx σσσσ ∑ − tti11 =1 i i ≤t este o rădăcină a lui σ(x) putem scrie: Pentru că Xi , 1≤ 1− tt ⋅++⋅+XXXσσσ, i = 1÷t (5 82) 0 = 11+−titii k i sumându-le obţinem ecuaţia: Înmulţind ecuaţiile (5 82) pe rând cu XiYi ş ttt kktkt 11+−++ ∑∑+⋅⋅++⋅+⋅XYXYXYσσ iitiiii 11∑− iii 111=== (5 83) t k 0∑=⋅⋅+XYσ iit i 1= sau, ţinând cont de (5 80) pentru k luând valorile 1,2, ,t: St+k + σ St+k-1 +…+ σ⋅Sk-1 +σ⋅Sk =0 , k = 1÷t (5 84) 1⋅t-1t Ecuaţiile (5 84) reprezintă un sistem de t ecuaţii cu t necunoscute (coeficienţii σ i) a cărui rezolvare constituie obiectivul acestui pas algoritmic Ecuaţiile (5 84) pot fi puse sub forma compactă: As⋅σ = Bs (5 85) unde: σS    1+t1 −ttSSS 11     2+tS2 +ttSSS 21σ     ; = σB (5 86) = s sA; =          tS2 −−tttSSS 2212tσ     Calculând inversa matricii As gă sim soluţia sistemului (5 85) în forma: -1 = (5 87) Aσ B ⋅ss Obs: Toate calculele trebuiesc făcute în câmpul GF(2²), atât coeficienţii sindrom, Sj, cât şi coeficienţii σ fiind elemente ale respectivului câmp În rezolvarea ecuaţiei (5 85) pot apărea trei situaţii: 1° rangul matricii As este e 0 5 u(i)=1; else u(i)=0; end end %Interliverul ina=1:N; for i=1:N y=(N-i+1)*rand(1); y=floor(y)+1; inb(i)=ina(y); if i N-3 & s0==s1 % readucere la zero u(i+j*N-N)=0; elseif i>N-3 u(i+j*N-N)=1; end v(1,i+j*N-N)=u(i+j*N-N); % primul bit din v af1=u(i+j*N-N)+s0+s1; % reactia codorului 1 af2=floor(af1/2); af1=af1-2*af2; if af1==s1 % iesirea codorului 1 af2=0; else af2=1; end s1=s0; s0=af1; v(2,i+j*N-N)=af2; % bitul livrat (via I) codorului 2 end end %Codorul 2 for j=1:n for i=1:N w(rnum(i)+j*N-N)=v(2,j*N-N+i); % interlivare end end for j=1:n s0=0; s1=0; % starea initiala for i=1:N v(2,i+j*N-N)=w(i+j*N-N); % bitul al doilea din v af1=w(i+j*N-N)+s0+s1; % reactia codorului 2 af2=floor(af1/2); af1=af1-2*af2; if af1==s1 % iesirea codorului 2 af2=0; else af2=1; end s1=s0; s0=af1; v(3,i+j*N-N)=af2; % bitul al treilea din v A-19 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta end end %Zgomotul B=10^(0 1*DB); sigma=1/(sqrt(2*R*B)); % dispersia zgomotului sig2=2*sigma*sigma; % dublul puterii zg %Semnalul receptionat for i=1:n*N for j=1:3 if v(j,i)==0 v(j,i)=-1; % modularea lui "0" prin "-1" end y=sigma*randn(1); % esantion de zgomot cu dispersie=sig2 r(j,i)=v(j,i)+y; end end % Erori inainte de decodare ercan=zeros(1,n); for j=1:n for i=1:N if r(1,j*N-N+i) 0 & v(1,j*N-N+i)==-1 ercan(j)=ercan(j)+1; end end end ercan %DECODAREA for j=1:n % contorul blocului in decodare % probabilitatea starilor initiale pentru dec1 si dec2 P(1,1)=1; P(2,1)=0; P(3,1)=0; P(4,1)=0; %Valoarea initiala a informatiei extrinseci for i=1:N Pw20(i)=0 5; end %j % Initializarea buclei iteratiilor pt blocul j for t=1:iter A-20 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta %t %Pw20 % Decodorul 2 % Recurenta inainte (dec 2) for i=1:N pw20=exp(-(r(2,i+j*N-N)+1)*(r(2,i+j*N-N)+1)/sig2); pw21=exp(-(r(2,i+j*N-N)-1)*(r(2,i+j*N-N)-1)/sig2); pv30=exp(-(r(3,i+j*N-N)+1)*(r(3,i+j*N-N)+1)/sig2); pv31=exp(-(r(3,i+j*N-N)-1)*(r(3,i+j*N-N)-1)/sig2); pw20=pw20/(pw20+pw21); pw21=pw21/(pw20+pw21); pv30=pv30/(pv30+pv31); pv31=pv31/(pv30+pv31); Px20(i)=pw20; P(1,i+1)=Pw20(i)*pw20*pv30*P(1,i)+(1-Pw20(i))*pw21*pv31*P(2,i); P(2,i+1)=(1-Pw20(i))*pw21*pv30*P(3,i)+Pw20(i)*pw20*pv31*P(4,i); P(3,i+1)=Pw20(i)*pw20*pv30*P(2,i)+(1-Pw20(i))*pw21*pv31*P(1,i); P(4,i+1)=(1-Pw20(i))*pw21*pv30*P(4,i)+Pw20(i)*pw20*pv31*P(3,i); sum=P(1,i+1)+P(2,i+1)+P(3,i+1)+P(4,i+1); P(1,i+1)=P(1,i+1)/sum; P(2,i+1)=P(2,i+1)/sum; P(3,i+1)=P(3,i+1)/sum; P(4,i+1)=P(4,i+1)/sum; end %P % Recurenta inapoi (dec 2) Pb(1,N+1)=P(1,N+1); Pb(2,N+1)=P(2,N+1); Pb(3,N+1)=P(3,N+1); Pb(4,N+1)=P(4,N+1); for i=N:-1:1 pw20=exp(-(r(2,i+j*N-N)+1)*(r(2,i+j*N-N)+1)/sig2); pw21=exp(-(r(2,i+j*N-N)-1)*(r(2,i+j*N-N)-1)/sig2); pv30=exp(-(r(3,i+j*N-N)+1)*(r(3,i+j*N-N)+1)/sig2); pv31=exp(-(r(3,i+j*N-N)-1)*(r(3,i+j*N-N)-1)/sig2); pw20=pw20/(pw20+pw21); pw21=pw21/(pw20+pw21); pv30=pv30/(pv30+pv31); pv31=pv31/(pv30+pv31); Pb(1,i)=Pw20(i)*pw20*pv30*Pb(1,i+1)+(1-Pw20(i))*pw21*pv31*Pb(3,i+1); Pb(2,i)=Pw20(i)*pw20*pv30*Pb(3,i+1)+(1-Pw20(i))*pw21*pv31*Pb(1,i+1); Pb(3,i)=Pw20(i)*pw20*pv31*Pb(4,i+1)+(1-Pw20(i))*pw21*pv30*Pb(2,i+1); Pb(4,i)=Pw20(i)*pw20*pv31*Pb(2,i+1)+(1-Pw20(i))*pw21*pv30*Pb(4,i+1); A-21 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta sum=Pb(1,i)+Pb(2,i)+Pb(3,i)+Pb(4,i); Pb(1,i)=Pb(1,i)/sum; Pb(2,i)=Pb(2,i)/sum; Pb(3,i)=Pb(3,i)/sum; Pb(4,i)=Pb(4,i)/sum; Pw20(i)=P(1,i)*Pb(1,i+1)*pv30+P(2,i)*Pb(3,i+1)*pv30+P(3,i)*Pb(4,i+1)*(1- pv30)+P(4,i)*Pb(2,i+1)*(1-pv30); Pw21(i)=P(1,i)*Pb(3,i+1)*(1-pv30)+P(2,i)*Pb(1,i+1)*(1- pv30)+P(3,i)*Pb(2,i+1)*pv30+P(4,i)*Pb(4,i+1)*pv30; Pw20(i)=Pw20(i)/(Pw20(i)+Pw21(i)); if Pw20(i)>1-ELIM Pw20(i)=1-ELIM; elseif Pw20(i) 1-ELIM Pv20(i)=1-ELIM; elseif Pv20(i)<ELIM Pv20(i)=ELIM; A-23 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta end end %Pb %Pu0 %uh % Inf extinsec dec 1 interlivata for i=1:N Pw20(rnum(i))=Pv20(i); % inf de la dec1 catre dec2 interlivata end % Erori la prezentul pas de iterare nrer(t,j)=0; for i=1:N if v(1,i+j*N-N)~=uh(i+j*N-N) nrer(t,j)=nrer(t,j)+1; end end nrer % Pasul iterativ urmator end % Blocul urmator end %Date finale %Index utilixat % i, j uzuali % t contorul iteratiilor %Variabile utilizate % af1 parametru temporar in codoare % af2 parametru temporar in codoare % DB RSZ in dB % ELIM prag pt prob minima % iter numarul de iteratii % k lungimea de constrangere % M numarul de stari posibile % n numarul de blocuri de informatie % N lungimea unui bloc de informatie % R rata de codare % sigma dispersia zgomotului % sig2 dublul puterii zg % sum suma de probabilitati % s0 starea primei celulei (pt ambele codoare) % s1 starea celulei a doua (pt ambele codoare) %Vectori utilizati A-24 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta % rnum(N) functia interliver % tu(n*N) timpul la informatie % tv(n*N/R) timpul in canal % u(n*N) secventa de date (informatia) % v(3,n*N) secventa codata / iesirea codorului 1 % w(n*N) iesirea codorului 1 interlivata A-25 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta Anexa F Decodarea în frecvenţă a codurilor Reed-Solomon Algoritmul Berlekamp- Massey Descrierea teoretică a algoritmului decodării 33 Ipoteze: - câmpul GF(2) generat prin p(x) = x + x + 1 (1) - codul RS (7, 3, 2) - polinomul generator dintre: 23 A gA(x) = (x + 1)(x + α )(x + α)(x + α) = 4235256 = x + αx + αx + αx + α = (1 3 6 6 7) (2) 234 B gB(x) = (x + α )(x + α)(x + α)(x + α) = 43323 = x + αx + x + αx + α = (1 4 1 2 4) (3) - cuvântul emis v(x) = 0 (4) 432 - cuvântul eroare ε(x) = αx + αx = 1-n i ==xε(0 0 0 5 0 3 0) ∑ i =0i (5) - cuvântul recepţionat w(x) = v(x) + ε(x) = ε(x) (6) TFD Consideraţii teoretice: - Transformata Fourier Discretă →←ε E ki 1-n kik E===αεαε 0,k 1-n () ik∑ 0i= (7) 1-n 1 i-ik- EE===ααε i 0, 1-n () ki∑ n 0 =k (8) - polinomul erorii Λ(x) = (xXi1 + 1) (xXi2 + 1) = 2 = 1 + Λ 1x + Λ2x (9) - locatorii erorilor Xi1, Xi2: -1-1i1i2 (Xi2) Xi1 = α, Xi2 = α Λ(Xi1) = 0 = Λ (10) n - ecuaţia erorilor: Λ(x)E(x) = Γ(x)(x + 1) (11) datorită gradelor polinoamelor: Γ(x) = ax + b (12) -j Obs: (11) se datorează: α , cu j = i1, i2, e rădăcină a lui Λ(x), ec (10) A-26 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta -j α, cu j ≠ i1, i2, e rădăcină a lui E(x), ec (8)(ej ≠ 0 doar pentru j = i1 şi i2) Soluţionarea problemei: - se identifică coeficienţii ecuaţiei (11): 0 x E0 = b 1 x E1 + Λ 1E0 = a 2 x E2 + Λ 1E1 + Λ2E0 = 0 3 x E3 + Λ 1E2 + Λ2E1 = 0 4 x E4 + Λ 1E3 + Λ2E2 = 0 5 x E5 + Λ 1E4 + Λ2E3 = 0 6 x E6 + Λ 1E5 + Λ2E4 = 0 7 x Λ 1E6 + Λ2E5 = b 8 x Λ 2E6 = a (13) - se află coeficienţii sindromului: 240 ()αααα =+= wS =  0 371   ααα () =+= wSα =  1   6432 () ⇒  B wAα αα =+= = Sα 2   563  ααα =+= wSα () =  3  624  = =+=ααα 1wS4 ()  = SEα =  00  E =  = α SE1 =  11  6  = B EA2α  ⇒ = SE = 22   E =  = α SE3 =  33   E1S = = 44  (14) i Λ din ec (13) - se află coeficienţii Λ 1 ş2 2  x E2 + Λ 1E1 + Λ2E0 = 0 3 A x E3 + Λ B 1E2 + Λ2E1 = 0  4 x E4 + Λ 1E3 + Λ2E2 = 0  + +   41323021EEEEEEEE ==Λ ==Λ 11 2121   + +   312201EEEEEE A B   22 + + 43424231EEEEEE   = =Λ = =Λ αα 2222   + + 312201EEEEEE   (15) - se află poziţiile eronate: A-27 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta 42 Λ(x) = α x + x + 1 (16) cu metoda descompunerii în factori a lui Λ(x): 4233 Λ(x) = α x + (α + α)x + 1 = α(αx + 1)x + αx + 1 = 3 = (α x + 1)(αx + 1) (17) -143 ⇒ xi1 = α ⇒ i1 = 3 ⇒ xi1 = α -16 xi2 = α ⇒ xi2 = α ⇒ i2 = 1 (18) cu metoda căutării: 044 ) = α + 1 + 1 = α ≠ 0 Λ(α 64 Λ(α) = α + α + 1 = α ≠ 0 225 Λ(α ) = α + α + 1 = α ≠ 0 333 Λ(α ) = α + α + 1 = 1 ≠ 0 454 Λ(α ) = α + α + 1 = 0 ⇒ poziţie eronată i1 = 7 – 4 = 3 555 Λ(α ) = 1 + α + 1 = α ≠ 0 626 Λ(α ) = α + α + 1 = 0 ⇒ poziţie eronată i2 = 7 – 6 = 1 (19) i v1 - se află valoarea caracterelor eronate v3 ş Se calculează polinomul Γ(x): A b = E0 = α a = E1 + Λ + α = 0 ⇒ Γ(x) = α 1E0 = α (20) 64 + α)/α = α B b = E0 =1/Λ 2(E2 + Λ1E1) = (α a = E1 + Λ Γ(x) = α 1E0 = E1 + Λ1b = 0 ⇒ (21) 3 x + 1)(αx + 1): Se derivează formal polinomul Λ(x) = (α ’333 Λ (x) = α(αx + 1) + α(αx + 1) = α+ α = 1 (22) Se calculează coeficienţii cuvântului eroare cu formula: -ii () ααΓ ε i i = i2 = pentru i = i1 ş i-i () 'αΛ (23) Obs: formula (23) se află derivând formal (11) şi ţinând cont de (8) şi de faptul că -i1 -i2 Λ(α) = 0 = Λ(α): ’’’nn-1 (x)E(x) + Λ (x)E(x) = Γ(x)(x + 1) + nxΓ(x) Λ -i -i -in cu x = α , i = i1 sau i = i2 ⇒ Λ(α) = 0 = (α) + 1 iar (24)devine: ’-i-i-nii-i Λ (α)E(α) = nαα Γ(α) de unde rezultă (23) A-28 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta 3444 ε α /1= α ⇒ v3 = w3 + ε = α + α = 0 3 = α3 222 ε /1= α v1 = w1 + ε = α + α = 0 1 = αα1 (24) Algoritmul Berlekamp-Massey –aplicaţie Cazul A : g(x) = gA(x) - se calculează coeficienţii sindrom: E0 = S0 = α E1 = S1 = α 6 E2 = S2 = α E3 = S3 = α - se află Λ(x) şi Γ(x) prin aplicarea algoritmului B-M: (0)(0) iniţializare r = 0 Λ(x) = 1 Γ(x) = 0 L0 = 0 (0)(0)-1 B(x) = 1 A(x) = x pasul r = 1 - se calculează discrepanţa: 6 E1 = α E2 = α E3 = α E(x) : E0 = α (0)000 Λ (x) = 1: Λ = 0 Λ = 0 Λ = 1 210 ∆ = 1 ⋅ E0 = α ≠ 0 - se calculează coef  1 dacă ∆ 0 şi 2L0 ≤ 0 = r - 1 1 ≠ δ 1   0 în rest ⇒ δ 1 = 1 -se calculează matricea D1: ∆ x1xE1x1α       01 = = = D       − 61-11 () −∆ δ 00Ex1α 011δ       ’’ (x) şi B(x): -se calculează Λ 0 + () + Λ  ()   Λ     0x1xE1xx'α = = =        61-01D () () 0ExBx'Bα        ’’ (x) şi A(x): -se calculează Γ 0 () ΓExx'α  () Γ    0 = == D     01 () 00xAx'A ()     -se calculează coeficientul L1: L1 = max(L0, r – L0) = max(0, 1) = 1 pasul r = 2 - se calculează discrepanţa: 6 E1 = α E2 = α E(x): E0 = α E3 = α ’’’’ Λ (x) = 1 + αx = 1 + E0x: Λ = 0 Λ = α Λ = 1 210 ’’224 ∆ + α = α ≠ 0 2 = Λ0E1 + Λ1E0 = E1 + E0 = α 0 şi 2L1 ≤ 1 - se calculează coef  1 dacă ∆ 2 ≠ A-29 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta δ 2   0 în rest ⇒ δ 2 = 0 -se calculează matricea D2: 42 ∆  () + x1xEE1x1α  2 01 = == D   12 − () −∆ x0x0x1δ 222δ   (2)(2) (x) şi B(x): -se calculează Λ -112 ()()   x1xEE1xx +   + ()() Λ Λ  10 = D = =      61-122 ()() xxExBxBα ()()    0   (2)(2) (x) şi A(x): -se calculează Γ 12 ()() Exxα  ()() ΓΓ  0 = ==D   ()() 122 00xBxA ()()   -se calculează coeficientul L2: L2 = max(L1, r – L1) = max(1, 1) = 1 pasul r = 3 - se calculează discrepanţa: 6 E1 = α E2 = α E3 = α E(x) : E0 = α (2)-1 Λ (x) = 1 + x = 1 + E0E1x: (2)(2)-1(2) Λ 2 = 0 Λ1 = E0E1 = 1 Λ0 = 1 (2)(2)(2)-12 65 ∆ + α = α ≠ 0 3 = Λ2E0 + Λ1E1 + Λ0E2 = E2 + E0E1= α - se calculează coef  1 dacă ∆ 0 şi 2L2 = 2 ≤ r – 1 = 2 3 ≠ δ 3   0 în rest ⇒ δ 3 = 1 -se calculează matricea D3: 52-1 ∆   () +   x1xEEE1x1α 3 102 = = = D − 13     − 221-1 () −∆ () + 00EEEx1α 333δ    102δ  (3)(3) (x) şi B(x): -se calculează Λ 22-1-1-1-13 ()     xEEEExEE1xEE1x () +++ + () Λ 10201010 D = = = 13 −     21-1-1-3 () ()( ) EEExEE1xExB () ++   102100   24   xx1α ++ =   22 xαα +   (3)(3) -se calculează Γ (x) şi A(x): 3 ()   () Γ   00EExα = = = −     3121-33D () () + ()   1020EEEE0xBα   -se calculează coeficientul L3: L3 = max(L2, 3 – L2) = max(1, 2) = 2 pasul r = 4 - se calculează discrepanţa: A-30 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta 6 E(x): E0 = α E1 = α E2 = α E3 = α (3)-1-22242 Λ (x) = 1 + E0E1x + E0(E0E2 + E1)x = 1 + x + αx : (3)-22(3)-1(3) Λ 2 = E0(E0E2 + E1) Λ1 = E0E1 Λ0 = 1 -1-1-23-23 ∆ 4 = E3 + E0E1E2 + E0E1E2 + E0E1 = E3 + E0E1= 0 - se calculează coef δ 4 = 0 - se calculează matricea D4: 3-2  () xEEE1 + 103 = D  4 x0  (4)(4) (x) şi B(x): -se calculează Λ 22-2-14 ()    () +++ () Λ 120010xEEEExEE1x = = D −    1-1244 () ++ ()() ()   101200xEE1EEEExB  2   + 23123021EEEEEEE + ++ 24xx1   22   xx1EEEEEEα ++ + + 120120   = =   222 α +   210xxEEα   + xx 22   ++ 120120EEEEEE   (4)(4) (x) şi A(x): -se calculează Γ 32  + 130EEE xEE +   () 2004  () xα Γ  +  1202EEE  ED == =   3204 ()   () 02ExAα   xEEE +    2120 + 120EEE  Polinoamele Λ(x) şi Γ(x) sunt: 2 ++EEEEEEE 24231230214 () () α++=++=Λ=Λxx1xx1xx 32 ++EEEEEE 120120 32 +EEE 4 130 () () α=+=Γ=ΓxExx 30 +EEE 120 utare: - se află poziţiile eronate ij prin metoda de că -i Λ(α ) = 0 i = 0, n-1 i i2 = 1 rezultă (vezi ec (19)): i1 = 3 ş i v1 - se află valorile caracterelor eronate: v3 ş - derivata formală a polinomului Λ(x) este (vezi (22)): Λ’(x) = 1 - coef cuvântului eroare sunt (vezi (23) şi (24)): 42 ε ε = α 3 = α1 - valorile caracterelor eronate sunt: v3 = w3 + ε 3 = 0 v1 = w1 + ε 1 = 0 - se selectează caracterele de informaţie: - cuvântul de cod este v = 0 0 0 0 0 0 0 - cuvântul de informaţie este i = 0 0 0 A-31 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta Cazul B g(x) = gB(x) - se calculează coeficienţii sindrom: E1 = S1 = α 6 E2 = S2 = α E3 = S3 = α E4 = S4 = 1 - se află Λ(x) şi Γ(x) prin aplicarea algoritmului B-M: (0)(0) iniţializare r = 0 Λ(x) = 1 Γ(x) = 0 L0 = 0 (0)(0) B(x) = 1 A(x) = 1 pasul r = 1 - se calculează discrepanţa: 6 E2 = α E3 = α E(x): E1 = α E4 = 1 (0)(0)(0)(0) Λ (x) = 1: Λ = 0 Λ = 0 Λ = 1 210 (0) ∆ ≠ 0 1 = E1Λ0 = α - se calculează coef  1 dacă ∆ 0 şi 2L0 ≤ 0 = r - 1 1 ≠ δ 1   0 în rest ⇒ δ 1 = 1 -se calculează matricea D1: ∆ x1xE1x1α      1 = D = =      − 61-11 () δ −∆ 00Ex1α  111δ     (1)(1) (x) şi B(x): -se calculează Λ () 1 + +  () Λ     1x1xE11xα = ==      61-11D () ()   1E1xBα    (1)(1) (x) şi A(x): -se calculează Γ () 1 xxE0xα   () Γ  1 = == D   1  () 1 001xA ()    -se calculează coeficientul L1: L1 = max(L0, r – L0) = max(0, 1) = 1 pasul r = 2 - se calculează discrepanţa ∆ 2: 6 E(x): E1 = α E2 = α E3 = α E4 = 1 (1)(1)(1)(1) (x) = 1 + E1x : Λ = 0 Λ = E1 Λ = 1 Λ 210 (2)(2)2 ∆ 0 2 = E1Λ1 + E2Λ0 = E1 + E2 = 1 ≠ - se calculează coef  1 dacă ∆ 0 şi 2L1 ≤ r – 1 = 1 2 ≠ δ 2   0 în rest ⇒ δ 2 = 0 -se calculează matricea D2: A-32 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta 2 ∆ x1xEE1x1   () +  212 = D = =     − 12 () −∆ x0x0x1δ  222δ    (2)(2) (x) şi B(x): -se calculează Λ () 5-12 +      + () + Λx1xEE1xE1xα  1 21 = D ==       ()61-1-22 () xxEExBα 1     1   (2)(2) -se calculează Γ (x) şi A(x): () 2 xxExα   () Γ  1 = = D   2  () 2 00xA ()    -se calculează coeficientul L2: L2 = max(L1, r – L1) = max(1, 1) = 1 pasul r = 3 - se calculează discrepanţa ∆ 3: 6 E(x): E1 = α E2 = α E3 = α E4 = 1 (2)-1(2)(2)-1E(2) Λ (x) = 1 + E1E2x: Λ = 0 Λ = E12 = 1 Λ = 1 210 -12 2 ∆ ≠ 0 3 = E3 + E1E2= α - se calculează coef  1 dacă ∆ 0 şi 2L2 ≤ r – 1 = 2 3 ≠ δ 3   0 în rest ⇒ δ 3 = 1 -se calculează matricea D3: 22-1  ∆   () +   x1xEEE1x1α 3 213 = D = = 13 −      − 521-1 () −∆ () + 00EEEx1α 333δ     213δ  (3)(3) (x) şi B(x): -se calculează Λ 22-1-1-1-13 ()     xEEEExEE1xEE1x () +++ () + Λ 21312121 D = = = −− 113     21-1-21-1-3 () ()() xEEEEEEEExExB () +++  213212131    25   xx1αα ++ =   35 xαα +   (3)(3) -se calculează Γ (x) şi A(x): xE   13 () xxExα  () Γ    21   xED = = =      6133 ()   x0xAα () 2      +EEE   213   -se calculează coeficientul L3: L3 = max(L2, r – L2) = max(1, 2) = 2 pasul r = 4 - se calculează discrepanţa ∆ 4: A-33 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta 6 E2 = α E3 = α E4 = E(x): E1 = α 1 (3)-1-1-122 Λ (x) = 1 + E1E2x + E1(E3 + E1E2)x: (3)-1-12(3)-1(3) Λ 2 = E1(E3 + E1E2) Λ1 = E1E2 Λ0 = 1 -1-1-126 ∆ ≠ 0 4 = E4 + E1E2E3 + E1E2(E3 + E1E2) = α - se calculează coef  1 dacă ∆ 0 şi 2L3 = 4 ≤ r – 1 = 3 4 ≠ δ 4   0 în rest ⇒ δ 4 = 0 - se calculează matricea D4: 63-2 ∆    () x1xEEE1x1α + 4 214 = D ==     − 14 () −∆ x0x0x1δ 444δ     (4)(4) (x) şi B(x): -se calculează Λ 22-1-1-14 ()   () +++ () Λ 213121xEEEExEE1x == 114D −−   21-1-21-4 () ()() +++ () 21321213xEEEEEEEExB   2  + + 2424233241EEEEEEE xx1α + +  ++  22xx1 = = EEEEEE + +  235231231  α + −−  2121-1-121-xxα  () () +++ 21321213xEEEEExEEE  (4)(4) -se calculează Γ (x) şi A(x): 32   + EEE 241 + xExE     252114 ()   () α + ΓxxEEEx2α +     xED231 = = =   ()62144     () xExA12α   xEEE +     2231   + EEE 231   Polinoamele Λ(x) şi Γ(x) sunt: 244 () () xx1xxα++=Λ=Λ 254 () () xxxxαα+=Γ=Γ utare: - se află poziţiile eronate ij prin metoda că -i Λ(α ) = 0 i = 0, n-1 i i2 = 1 rezultă (vezi ec (19)): i1 = 3 ş i v1 - se află valorile caracterelor eronate: v3 ş - derivata formală a polinomului Λ(x) este (vezi (22)): Λ’(x) = 1 - coef cuvântului eroare sunt (vezi (23) şi (24)): 3-3-33-2-134 ε Γ(α) /Λ’(α) = α(α +α) = α⋅α = α 3 = α -1-132 ε (α) /Λ’(α) = α(1 +α) = α⋅α = α 1 = αΓ - valorile caracterelor eronate sunt: v3 = w3 + ε 3 = 0 A-34 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta v1 = w1 + ε 1 = 0 - se selectează caracterele de informaţie: - cuvântul de cod este v = 0 0 0 0 0 0 0 - cuvântul de informaţie este i = 0 0 0 A-35 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta Anexa G 1 Calculul funcţiei de transfer T(δ,β,λ), relaţia (6 3) Utilizând notaţiile din Figura G 1, se pot scrie ecuaţiile: 2  ⋅⋅+⋅⋅⋅=vλβxλβδw  ⋅⋅+⋅⋅⋅=uλβδwλβδu (G 1) ⋅   λδwλδv⋅⋅+⋅⋅=u  βλ δ u βλ δλ δ 22 λ βλ δ δ δλ x wvy βλ Figura G 1 Diagrama de stare a codorului convoluţional din Figura 6 1 Ultimele două ecuaţii din (G 1) se pot scrie: u = β⋅v δ⋅λ⋅w = (1 + δ⋅β⋅λ)⋅v Eliminând pe w între ultima ecuaţie şi prima din (G 1) rezultă: 232 v⋅(1 + δ⋅β⋅λ + δ⋅β⋅λ ) = δ⋅β⋅λ⋅x 2 Din Figura G 1 rezultă că y = δ ⋅λ⋅v, ca atare putem calcula funcţia de transfer: 35 λβδy⋅⋅ () λβ,δ,T== (G 2) () λ1δβλ1x+⋅+ Ecuaţia (G 2) se poate obţine şi prin regula lui Mason: A-36 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta −⋅)B(1C ∑∑ iji i T = (G 3) −B1 ∑ j j a căii „i” între nodurile final şi iniţial; unde: –Ci = transmitanţ –Bj = „transmitanţ a” buclei „j”, dintre cele ne-adiacente, peste tot graful; –Bij = „transmitanţ a” buclei „j”, dintre cele ne-adiacente căii „i”; Astfel, în graful din Figura G 1 există două bucle ne-adiacente, una aferentă nodului „u”, de transmitanţă B1 = δ⋅β⋅λ , iar cea de-a doua între nodurile „w” şi „v”, de transmitanţă B2 2 De asemenea există două căi între nodurile „x” şi „y”, şi anume: x-w-u-v-y şi x- = δ⋅β⋅λ 62453 ⋅β⋅λ şi C2 = δ⋅β⋅λ Cu acestea, funcţia de transfer w-v-y, având transmitanţele C1 = δ se scrie: 42635 +⋅+ )B(1CC⋅ βδλ)βδ(1λβδ⋅+⋅⋅+⋅⋅⋅ 221λ T(δ,β,λ) = = 2 ++ βδλβδ1BB1⋅ ⋅+⋅⋅+ 21λ 35 λβδ⋅⋅ = () λ1δβλ1+⋅+ Trebuie menţionat că, în cazul sumei modulo-doi, operaţiile de adunare şi scădere sunt identice Vectorul recepţionat: 01 00 01 00 00 0 00 00000000 11 11111111 111111 000000 10101010 Tactul 0 010101 01 010101 101010 Figura G 2 Diagrama trellis şi semnalul recepţionat în cazul canalului binar simetric A-37 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta 2 Exemplificarea algoritmului Viterbi În cazul canalului CBS, diagrama trellis şi semnalul recepţionat, arată ca în Figura G 2 Algoritmul recepţiei presupune, la fiecare tact, calculul metricii (relaţia (6 14)) şi selecţia căii (ramurii) celei mai plauzibile Deoarece în primele două tacturi, diagrama trellis nu conţine noduri cu două intrări, selecţia se face începând cu tactul trei, Figura G 3a Vectorul recepţ 01 00 01ionat: 0 –1 –1 –2 00 0000 –3 11 1111 –1 –3 –2 11 00 –3 10 10–2 –5 Tactul 3 01 01 –2 01 –2 –3 a) 10 –4 Vectorul recepţionat: 01 00 01 –2 00 0000 11 1111 –2 Tactul 3 01 01 –2 01 b) –3 Figura G 3 Diagrama trellis până la tactul 3, a) –înainte şi b) –după selecţia aferentă tactului trei; După selecţia de la tactul 3, au rămas doar patru drumuri posibile, având distanţa faţă de secvenţa recepţionată –2, (pentru primele trei), şi –3 (pentru ultimul) Aceste A-38 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta distanţe constituie „zestrea” fiecărei căi şi, totodată, datele pentru un nou pas algoritmic, cel aferent tactului 4 Figurile G 4 şi G 5 prezintă operaţiile aferente tacturilor 4 şi 5 Vectorul recepţ ionat: 01 00 01 00 –2 –2 00 000000 –4 11 111111 –2 –4 11 00–2 10 –3 Tactul 4 0101 01 –2–4 0101 a) –3 –3 10 –4 Vectorul ionat: recepţ01 00 01 00 –2 00 000000 11 11 00 –2 10–3 Tactul 4 01 01 b) 01 –3 Figura G 4 Diagrama trellis la tactul 4: a) –înainte şi b) –după selecţie; Vectorul recepţionat: 01 00 01 00 00 –2 –2 00 00000000 A-39–5 11 1111 –4 11 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta a) Vectorul ionat: recepţ01 00 01 00 00 –2 00 00000000 11 11 0000 –3 1010 –3 Tactul 5 01 01 b) 01 –3 Figura G 5 Diagrama trellis la tactul 5: a) –înainte şi b) –după selecţie; Acelaşi algoritm, exemplificat pentru canalul AWGN, este prezentat prin figurile următoare Metrica ramurii se calculează în acest caz după relaţia (6 16), iar ramura selectată este cea cu metrică mai mare (dintre cele două ce converg în acelaşi nod) Vectorul ionat: recepţ+1,+5 –1,+3 +2, –3 +3,+6 +7,–1 0 +1,+1 +1,+1+1,+1+1,+1+1,+1 A-40 -1,-1-1,-1-1,-1 -1,-1 -1,-1 -1,-1-1,-1-1,-1 1111 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta a) Vectorul ionat: recepţ+1,+5 –1,+3 +2, –3 0 6 8 7 +1,+1 +1,+1 +1,+1 –1 -1,-1-1,-1-1,-1 –6 4 9 -1,-1 +1,+1 –3 b) înainte de selecţie -1,+1–2 –1 -1,+1 Tactul 3 +1,-1+1,-1 –5 +1,-1 –10 9 Vectorul -1,+1 recepţionat: –1,+3 +2, –3 –15+1,+5 7 +1,+1 +1,+1 +1,+1 -1,-1-1,-1 9 c) după selecţie -1,+1–1 Tactul 3 +1,-1 9 Figura G 6 Diagrama trellis în cazul canalului AWGN A-41 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta Vectorul recepţionat: +1,+5 –1,+3 +2, –3 +3,+6 7 16 +1,+1 +1,+1+1,+1+1,+1 –10 -1,-1 -1,-1-1,-1 9 –2 -1,-1 +1,+1 8 d) înainte de selecţie -1,+1-1,+1–1 12 Tactul 4 +1,-1 +1,-1 6 +1,-1 9 6 -1,+1 12 Vectorul ionat: +1,+5 –1,+3 +2, –3 +3,+6recepţ 16 +1,+1 +1,+1+1,+1+1,+1 -1,-1 -1,-1 +1,+1 8 e) după selecţie -1,+1-1,+1 12 Tactul 4 +1,-1 -1,+1 12 Figura G 6 Diagrama trellis în cazul canalului AWGN A-42 Coduri utilizate în sistemele de transmisiuni cu spectru împrăştiat Referat nr 1 Horia Balta Vectorul ionat: recepţ+1,+5 –1,+3 +2, –3 +3,+6 +7,–1 16 22 +1,+1 +1,+1+1,+1+1,+1+1,+1 6 -1,-1 -1,-1-1,-1 10 -1,-1 +1,+1+1,+1 8 18 f) înainte de selecţie -1,+1-1,+1-1,+112 0 Tactul 5 +1,-1+1,-1 20 +1,-1 16 -1,+1-1,+1 12 4 Vectorul ionat: recepţ+1,+5 –1,+3 +2, –3 +3,+6 +7,–1 22 +1,+1 +1,+1+1,+1+1,+1+1,+1 -1,-1 -1,-1 +1,+1+1,+1 18 g) după selecţie -1,+1-1,+1 Tactul 5 +1,-1+1,-1 20 +1,-1 16 -1,+1 Figura G 6 Diagrama trellis în cazul canalului AWGN A-43